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THEORETICAL  STUDY  OF  THE  TRANSONIC  LIFT  OF  A  DOUBLE -WEDGE ' 
PROFILE  WITH  DETACHED  BOW  WAVE  1 
By  Walter  G.  Vincent!  and  Cleo  B.  Wagoner 


SUMMARY 


A  theoretical  study  is  described  of  the  aerodynamic  characteristics 
at  small  angle  of  attack  of  a  thin,  double -wedge  profile  in  the  range 
of  supersonic  flight  speed  in  which  the  bow  wave  is  detached.  The 
analysis  is  carried  out  within  the  framework  of  the  transonic  (nonlinear) 
small -disturbance  theory,  and  the  effects  of  angle  of  attack  are  regarded 
as  a  small  perturbation  on  the  flow  previously  calculated  at  zero  angle. 
The  mixed  flow  about  the  front  half  of  the  profile  is  calculated  by 
relaxation  solution  of  a  suitably  defined  boundary-value  problem  for  the 
transonic  small -disturbance  equation  in  the  hodograph  plane  (i.e.,  the 
Tricomi  equation) .  The  purely  supersonic  flow  about  the  rear  half  is 
found  by  an  extension  of  the  usual  numerical  method  of  characteristics. 
Analytical  results  are  also  obtained,  within  the  framework  of  the  same 
theory,  for  the  range  of  speed  in  which  the  bow  wave  is  attached  and  the 
flow  is  completely  supersonic. 

The  calculations  provide,  for  vanishingly  small  angle  of  attack, 
the  following  information  as  a  function  of  the  transonic  similarity 
parameter:  (l)  chordwise  lift  distribution,  (2)  lift-curve  slope,  and 
(3)  position  of  center  of  lift.  As  in  previous  studies,  the  aerodynamic 
characteristics  of  a  profile  of  given  thickness  ratio  show  little  varia¬ 
tion  with  free -stream  Mach  number  as  the  Mach  number  passes  through  1. 

As  the  Mach  number  is  increased  to  higher  values,  however,  the  lift- 
curve  slope  rises  to  a  pronounced  maximum  in  the  vicinity  of  shock 
attachment  and  then  declines.  Correspondingly,  the  center  of  lift  moves 
forward  toward  the  leading  edge  and  then  returns  aft.  These  findings 
are  in  marked  contrast  to  the  behavior  of  the  drag  coefficient  at  zero 
angle  of  attack,  which  was  found  in  earlier  work  to  decrease  monotoni- 
cally  as  the  Mach  number  increased  above  1.  At  Mach  numbers  above  that 
for  shock  attachment,  the  results  of  the  present  calculations  tend  toward 
those  given  by  classical  linear  theory. _ _  ■  _ • 

Portions  of  this  work  were  reported  at  the  VUIth  International 
Congress  on  Theoretical  and  Applied  Mechanics,  Istanbul,  Turkey, 

August  20-28,  1952. _ 
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INTRODUCTION 

The  theoretical  problem  of  the  transonic  flow  over  a  thin,  double -' 
wedge  profile  at  zero  angle  of  attack  has  been  treated  in  several  papers 
in  recent  years.  These  papers  have  in  common  that  they  employ  the 
simplifying  concepts  of  the  transonic  small -disturbance  theory  and 
utilize  the  holograph  transformation  to  linearize  the  resulting  mathe¬ 
matical  problem.  Following  this  approach,  Guderley  and  Yoshihara 
(reference  l)  began  by  solving  the  problem  for  a  free -stream  Mach 
number  of  1,  using  analytical  methods  for  the  mixed  flow  over  the  front 
wedge  and  the  method  of  characteristics  for  the  purely  supersonic  flow 
over  the  rear.  Somewhat  later,  the  present  authors,  using  a  combina¬ 
tion  of  relaxation  methods  and  the  method  of  characteristics  (refer¬ 
ences  2  and  3)>  extended  the  results  to  free -stream  Mach  numbers  greater 
than  1,  where  a  detached  bow  wave  occurs  ahead  of  the  profile.  At  about 
the  same  time.  Cole  (reference  4)  obtained  an  analytical  solution  for 
the  flow  over  the  front  wedge  at  subsonic  flight  speeds,  utilizing,  in 
effect,  the  special  assumption  of  a  vertical  sonic  line  from  the 
shoulder  of  the  wedge.  More  recently.  Trilling  (reference  5)  has  been 
able  to  remove  this  special  assumption  and,  with  the  aid  of  less  strin¬ 
gent  approximations  regarding  the  flow  over  the  rear  wedge,  to  extend 
the  solution  for  the  subsonic  case  to  Include  the  complete  profile.  As 
a  result  of  these  investigations,  the  problem  of  the  double-wedge 
profile  at  zero  angle  of  attack  may  be  regarded  as  substantially  solved 
within  the  limitations  of  the  transonic  small-disturbance  theory.  The 
experimental  studies  of  Liepmann  and  Bryson  (references  6  and  7)  and 
Griffith  (reference  8)  indicate  that  the  theoretical  findings  are  in 
fundamental  agreement  with  the  physical  facts. 

In  a  recent  paper  (reference  9),  Guderley  and  Yoshihara  have  con¬ 
tinued  their  investigations  of  the  double -wedge  profile  at  Mach  number  1 
by  considering  the  influence  of  a  vanishingly  small  angle  of  attack. 

The  basic  idea  in  this  later  work  is  to  regard  the  effects  of  angle  of 
attack  as  a  first-order  perturbation  on  the  nonlinear  flow  previously 
calculated  at  zero  angle.  This  approach  leads  to  a  linear  boundary- 
value  problem  in  both  the  physical  and  hodograph  planes.  The  calcula¬ 
tion  for  the  front  wedge  is  still  carried  out,  however,  in  the  hodograph 
plane,  since  the  basic  procedures  can  then  be  taken  over  directly  from 
the  previous  work.  By  this  means,  Guderley  and  Yoshihara  obtain  results 
for  the  lift-curve  slope  of  the  profile  at  zero  angle  of  attack  and  for 
the  corresponding  distribution  of  lift  along  the  chord. 

The  aim  of  the  present  paper  is  to  extend  the  results  for  the 
double  wedge  at  angle  of  attack  to  the  case  of  supersonic  flight  with 
detached  bow  wave.  The  fundamental  ideas  of  Guderley  and  Yoshihara  are 
followed  in  reducing  the  calculations  for  the  front  wedge  to  a  perturba¬ 
tion  problem  in  the  hodograph  plane .  The  detailed  formulation  of  the 
problem  is,  however,  necessarily  different  in  the  present  case.  The 
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boundary  conditions  for  the  problem  appear  in  terms  of  the  results 
already  obtained  at  zero  lift  (references  2  and  3) ,  and  the  solution  is 
carried  out  by  numerical  methods  which  differ  only  slightly  from  those 
devised  for  the  earlier  work.  The  lift  on  the  rear  wedge  is  calculated 
by  an  extension  of  the  method  of  characteristics.  The  body  of  the 
paper  is  devoted  to  the  detailed  formulation  of  the  boundary- value 
problem  in  the  holograph  plane  and  to  a  discussion  of  the  final  results. 
Noteworthy  differences  between  the  numerical  procedures  used  in  the 
present  work  and  those  already  described  in  reference  3  are  treated  in 
appendices  at  the  end  of  the  report. 


NOTATION 
Primary  Symbols 


a*'  critical  speed  (i.e.,  speed  at  which  the  speed  of  flow  and 

the  speed  of  sound  are  equal) 

b  numerical  coefficient 

(See  equations  (39)  and  (4o).) 

c  airfoil  chord 

c7  lift  coefficient  ( l-1?1  Per  unit  sPan 

\  <loc 


-m 


moment  coefficient  for  moments  taken  about  leading 
edge  (  momen £  Per  unit  span  j 

6  86  \  qQc2  J 


pressure  coefficient 


(l  -  Po 

\  !o 


Iw  integral  defined  by  equation  (45) 

k 2  numerical  constant 

(See  equation  (10).) 


M  Mach  number 


m  slope  of  segment  of  Mach  line  in  characteristics  net 


P 


static  pressure 


4  NACA  TN  2832 

Ap  local  lifting  pressure  (i.e.,  difference  between  static 

pressures  on  upper  and  lower  surfaces) 

q  dynamic  pressure 

t  airfoil  thickness 

V  speed  of  flow  ‘ 

x,y  Cartesian  coordinates 


0 
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generalized  Cartesian  coordinates 
(See  equations  (43).) 

chordwise  position  of  center  of  lift 


slope  of  curve  of  lift  coefficient  versus  true  angle  of 
attack  evaluated  at  zero  angle 

slope  of  curve  of  moment  coefficient  versus  true  angle  of 
attack  evaluated  at  zero  angle 

normalized  angle  of  attack;  also  denotes  true  angle  of 

dC7 

attack  when  used  in  derivatives  such  as  - ,  etc. 

da  * 

absolute  value  of  rj  at  left-hand  limit  of  lattice 
ratio  of  specific  heats  (1.4  for  air)  , 


A  basic  lattice  interval 

t  function  of  q  and  9 

(See  equation  (A6) . ) 

•V 

Tbh  normalized  speed  of  flow 

(See  equation  (la)  and  page  10.) 

/  '' 

special  values  of  q 

(See  sketch  (m)  on  page  46-) 

«v» 

normalized  inclination  of  flow;  9  also  denotes  true  inclina¬ 
tion  of  flow  in  equation  (lb) 

(See  equation  (lb)  and  page  10.) 

8V  normalized  half  angle  of  wedge 

£0  transonic  similarity  parameter 

(See  equation  (13).) 
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p  fluid  density 

w  i|r  stream  function 

5\(f'A  incremental  values  of  stream  function 

(See  equations  (A 9)  and  (All).) 


Subscripts 


a,b,c  points  in  characteristics  net 

(See  page  page  43.) 

A,B  components  of  total  stream  function 

(See  equation  (39)-) 

o  conditions  in  free  stream 

s  singular  solution 

(See  equations  (a6)  and  (A7).) 

0,1,2,  value  at  a  prescribed  lattice  point 

etc. 

*  conditions  at  critical  speed 


Superscripts 


(  )  quantities  determined  at  zero  angle  of  attack 

(  ) *  derivative  with  respect  to  normalized  angle  of  attack 
evaluated  at  zero  angle 


* 
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BOUNDARY-VALUE  PROBLEM  IN  HODOGRAPH  PLANE 
Description  of  Flow  Field 

Sketch  (a)  is  a  drawing  of  the  idealized,  inviscid  flow  which 


Mach 

lines 


may  he  expected  about  a  wedge  profile  when  the  angle  of  attack  is 
sufficiently  less  than  the  semiapex  angle  of  the  wedge.  Sketch  (b) 
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shows  the  corresponding  hodograph  representation  of  the  flow  over  the 
front  wedge,  which  is  the  region  of  prime  theoretical  concern.  Except 
for  the  substitution  of  the  detached  how  wave  in  place  of  the  infinite 
free  stream,  these  representations  follow  the  lines  assumed  by  Guderley 
and  Yoshihara  in  reference  9-  The  corresponding  drawings  for  zero  angle 
of  attack,  which  are  fundamental  to  the  present  case,  have  been  described 
in  detail  in  references  2  and  3. 

In  the  present  example,  the  path  of  the  central  streamline  in  the 
physical  and  hodograph  planes  is  briefly  as  follows:  The  streamline 
leaves  the  bow  wave  in  the  physical  plane  (or  the  shock  polar  in  the 
hodograph  plane)  at  point  A.  It  then  proceeds  with  decreasing  subsonic 
speed  to  a  stagnation  point  0  on  the  underside  of  the  profile.  At  0 
the  streamline  branches.  The  lower  branch  runs  downstream  along  the 
lower  surface  of  the  profile  with  fixed  inclination  but  increasing 
speed.  The  sonic  speed  is  reached  at  the  shoulder  L,  where  the  speed 
then  increases  discontinuously  in  accord  with  the  Prandtl -Meyer  rela¬ 
tions.  The  shoulder  itself  maps  in  the  hodograph  onto  the  upgoing 
characteristic  EM.  The  upper  branch  of  the  central  streamline  proceeds 
from  0  upstream  along  the  surface  of  the  profile.  The  inclination  here 
is  again  fixed  by  that  of  the  surface,  and  the  speed  increases  to  the 
sonic  value  at  the  leading  edge  J.  At  this  point  the  flow  is  charac¬ 
terized  by  another  Prandtl -Meyer  expansion  to  supersonic  speed. 

The  flow  configuration  which  should  be  assumed  on  the  upper  surface 
near  the  leading  edge  is  open  to  conjecture.  Since  the  geometrically 
avaiiabie  angie  of  turn  will,  for  any  thin  airfoil,  be  greater  than  the 
130  permissible  for  expansion  to  a  vacuum,  a  region  of  separation  is 
to  be  expected.  If  the  angle  of  attack  is  not  too  great,  this  region 

Probably  be  closed,  with  the  central  streamline  reattaching  to  the 
upper  surface  a  small  distance  behind  the  leading  edge.  This  reattach¬ 
ment  will  be  followed  by  a  compression  of  the  flow  through  a  system  of 
shock  waves  whose  arrangement  is  sketched  only  formally  in  the  physical 
plane  (and  not  at  all  in  the  hodograph  plane,  where  the  correct  repre¬ 
sentation  would  probably  lie  on  several  sheets).  The  effects  of  the 
ilow  near  the  leading  edge  will  be  mentioned  later,  but  the  exact 
process  will  remain  undefined.  Whatever  the  details,  the  speed  on  the 
upper  surface  will  return  to  a  subsonic  value  at  some  point  K  just 
downstream  of  a  terminating,,  normal  shock  wave.  From  K  the  central 
streamline  continues  at  fixed  inclination  downstream  along'  the  upper 
surface,  the  speed  increasing  once  more  to  the  sonic  value  at  the 
shoulder  B.  At  this  point  another  expansion  takes  place,  similar  to 
that  which  occurs  at  the  corresponding  point  on  the  lower  surface.  In 
this  case  the  shoulder  is  represented  in  the  hodograph  by  the  down¬ 
going  characteristic  BG. 

The  supersonic  expansion  fan  from  the  shoulder  at  B  (and  similarly 
at  L)  is  discussed  in  detail  in  references  2  and  3.  Suffice  it  here  to 
say  that  the  supersonic  flow  field,  of  which  the  expansion  fan  is  the 
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initial  part,  is  separated  into  two  regions  by  the  Mach  line  GE,  which 
runs  from  the  shoulder  to  the  sonic  point  on  the  bow  wave.  (This  line 
was  termed  the  ’'separating"  Mach  line  in  reference  2.)  The  supersonic 
flow  in  the  region  upstream  of  the  Mach  line  GE  is  interdependent  with 
the  subsonic  field  between  the  bow  wave  and  sonic  line.  To  obtain  a 
solution  for  the  front  wedge,  a  problem  in  transonic  flow  must  there¬ 
fore  be  solved  for  the  subsonic  field  and  the  interdependent  portion  of 
the  supersonic  expansion  fan.  Conditions  in  the  supersonic  flow  down¬ 
stream  of  the  Mach  line  GE  have  no  influence  upon  the  subsonic  field. 
The  continuation  of  the  flow  beyond  GE  can  be  accomplished  by  purely 
supersonic  methods  once  the  solution  of  the  transonic  problem  is  known. 

Aside  from  the  obvious  lack  of  symmetry  in  the  present  case,  the 
main  difference  between  the  flow  here  and  that  previously  studied  at 
zero  angle  of  attack  is  the  existence  in  the  present  problem  of  the 
localized  supersonic  region  in  the  vicinity  of  the  nose.  As  pointed 
out,  conditions  in  this  region  are  difficult  to  formulate.  The  problem 
has  been  considered  by  Guderley  and  Yoshihara  (reference  9)  i*1  "the 
course  of  their  work  at  Mach  number  1.  They  find  that,  if  the  nose 
region  is  disregarded  in  the  hodograph  and  the  boundary  condition 
along  KB  is  fulfilled  all  the  way  into  0,  then  the  influence  on  the 
lift  of  the  resulting  fictitious  flow  at  the  nose  is  of  somewhat  higher 
than  the  second  order  in  the  angle  of  attack.  This  suggests  that  the 
effects  of  the  real  flow  at  the  nose  may  be  neglected  in  a  first-order 
analysis  such  as  the  present.  In  the  work  which  follows,  as  in  the 
calculations  of  Guderley  and  Yoshihara,  the  supersonic  region  at  the 
leading  edge  will  therefore  be  disregarded. 


Formulation  of  Boundary-Value  Problem 


As  in  reference  3,  the  analysis  is  based  on  the  equations  of  the 
transonic  small -disturbance  theory  with  the  stream  function  V  as  the 
dependent  variable.  If  the  effects  of  the  flow  at  the  nose  are  ignored, 
the  problem  of  the  wedge  at  angle  of  attack  a  is  then  readily  formu¬ 
lated  as  a  boundary-value  problem  in  the  hodograph  plane.  To  solve  this 
problem  for  vanishingly  small  jx,  it  will  be  assumed  that  the  solution_  'lf 
at  angle  of  attack  can  be  expressed  as  the  sum  of  the  basic  solution  V 
previously  obtained  at  zero  angle  plus  a  perturbation  term  a^f’ ,  where  i|r 
is  a  f -unction  which  does  not  itself  involve  a.  By  consideration  of  the 
difference  between  the  boundary-value  problems  for  f  and  V,  a  problem 
for  the  perturbation  function  can  be  formulated.  The  boundaries  for 

this  problem  turn  out  to  be  the  same  as  those  for  the  problem  at  zero 
angle,  and'  the  boundary  values  themselves  appear  in  terms  of  t .  The 
details  of  these  matters  will  now  be  given.  The  reader  who  is  interested 
only  in  the  results  can  turn  directly  to  the  section  Chordwise  Distribu¬ 
tion  of  Lift  on  page  28. 
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Basic  equations. -  The  basic  equations  will  be  taken  in  the  form 
given  in  reference  3>  that  is,  in  terms  of  small  disturbances  from  the  ^ 
critical  speed  a*.2  The  independent  variables  are  the  normalized  speed  T) 
arid  the  normalized  inclination  0  as  defined  by  the  relations 


V/a»-l 
~  V0/a* -  1 


(lb) 


where 

a 

V  local  speed  of  flow 

0  local  inclination  of  flow  relative  to  direction  of  free  stream 
V0  free -stream  speed 

a*  critical  speed  (i.e.,  speed  at  which  the  speed  of  flow  and  the 
speed  of  sound  are  equal) 

7  ratio  of  specific  heats 

Use  of  these  variables  is  equivalent  to  introducing  the  rules  for  tran¬ 
sonic  similarity.  In  terms  of  the  foregoing  holograph  variables,  the 
differential  equation  for  the  stream  function  \|r  as  given  by  the  tran¬ 
sonic  small-disturbance  theory  Is 


(2) 


2As  discussed  in  several  recent  papers  (e.g.,  references  10  and  ll),  the 
theory  can  also  be  formulated  in  terms  of  disturbances  from  the  free- 
stream  speed  VQ.  This  latter,  less  restrictive  formulation  reveals 
clearly  the  relationship  which  exists  between  the  transonic  small- 
disturbance  theory  and  the  familiar  linear  theory  of  subsonic  or  super¬ 
sonic  flow.  As  shown  by  Spreiter  (see  page  9  of  reference  ll),  an  a* 
analysis  will  yield  values  of  the  pressure  coefficient  identical  to 
those  of  a  VQ  analysis  provided  the  pressure  coefficient  and  similar¬ 
ity  parameter  in  the  former  case  are  taken  as  in  equations  (4)  and  (13) 
below.  If  this  procedure  is  followed,  the  results  of  the  a*  analysis 
may  even  be  expected  to  tend  toward  those  of  linear  theory  as  the  free- 
stream  Mach  number  increases  or  decreases  from  1.  (An  analytical 
example  of  just  this  behavior  has  been  given  by  Bryson  in  appendix  A 
of  reference  7-)  It  appears,  therefore,  that  the  a*  formulation, 
when  suitably  used,  gives  results  of  wider  theoretical  validity  than 
■would  be  anticipated  on  the  basis  of  its  own  rather  restrictive  under- 
lying  assumption. _ _ _  • 


♦tvh  "  2 ^  +: 


W  =  O 
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This  is  essentially  the  linear  differential  equation  first  studied  hy 
Tricomi  (reference  12).  It  is  elliptic  for  rf<0  (subsonic  speeds) 
and  hyperbolic  for  rf >0  (supersonic  speeds). 


The  transformation  from  the  hodograph  to  the  physical  plane  is 
governed  by  the  differential  relations 


dx  = 


1 

r(7+l)  (V0/a*-l) 

P*a* 

2 

1/2 


(2^?f  +  *^L~) 


(3a) 


dy  =  -i-  (WL~  +  V®)  =  7TT  dt 
P*a*  T*  “  o*a* 


(31) 


where  x=x(rf,0)  and  y=y(?f,§”)  are  physical  coordinates  (horizontal  and 
vertical,  respectively)  corresponding  to  a  given  velocity  r\}0.  The 
symbol  denotes  the  fluid  density  at  the  critical  speed  a*.  Within 

the  approximation  of  the  transonic  small -disturbance  theory,  the  pressure 
coefficient  Cp  =  (p-p0)/q0  can  be  calculated  from  the  relation 


Cp  =  -2  =  -2(V0/a* -l)(rf-l)  (*0 

a* 


The  local  Mach  number  is  related  to  the  speed  of  flow  by  the  equation 


M5-! 

7+1 


For  simplicity  of  notation,  the  tilde  will  be  omitted  from  the 
symbols  t]  and  0  in  the  remainder  of  the  development.  It  is  to  be 
understood,  unless  stated  otherwise,  that  the  quantities  r\  and  0  are 
themselves  the  normalized  quantities  defined  by  equations  (l). 


Problem  at  zero  angle  of  attack. -  When  the  angle  of  attack  is  zero, 
the  localized  region  of  supersonic  flow  at  the  leading  edge  disappears 
from  sketch  (a),  and  the  flow  field  becomes  symmetrical  about  the  chord 
line.  The  corresponding  boundary-value  problem  in  the  r \,Q'  plane  has 
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been  set  forth  in  reference  3.  It  is  restated  in  sketch  (c),  where  both 
the  upper  and  lower  halves  of  the  flow  field  are  now  included.  In  this 
representation,  the  surfaces  of  the  wedge  appear  as  the  semi -infinite 
horizontal  lines  OB  and  OL,  and  the  subsonic  portion  of  the  shock  polar 
appears  as  the  curve  ME. 


Sketch  (c) 

_  If  the  stream  function  foir  zero  angle  of  attack  is  denoted  by 

the  differential  equation  to  be  satisfied  here  is  given  by 
equation  (2)  as 

-  ^*ee>  -  0  («) 

The  requirement  that  the  flow  shall  be  tangent  to  the  surfaces  of  the 
wedge  provides  the  boundary  conditions 

^(■n#±Sw)  =0  for  T)  <  0  (7) 
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where  0W  denotes  the  normalized  half -angle  at  the  leading  edge.  The 
stagnation  point  at  the  leading  edge  is  represented  in  the  present 
theory  by  the  condition  that 

f — >0  for  T)— >-»,  -,&w'50  =  ew 

Along  the  shock  polar  ME,  the  relations  for  an  oblique  shock  wave 
require  that 

*  i  =0  (9) 

*1  3+5n  0 

for 

0  =  ±  (l-ii)  JT+r\  ,  -1  <  T)  <  0 


Along  the  sonic  line,  boundary  conditions  are  prescribed  which  represent 
the  influence  exerted  on  the  subsonic  field  by  the  interdependent  portion 
of  the  supersonic  expansion  fans.  On  the  basis  of  the  procedures  given 
in  reference  3,  this  influence  can  be  represented  completely  by  the 
requirement  that 


*9(0,%) 

[±(0!-0)]2/3 


(10) 


where  the  upper  signs  apply  for  1<0<0w  and  ttie  lover  siSns 
for  -0w  <  0  <  -1 .  The  constant  k2  which  appears  here  is  given  by 


_  24/3  if  ~ 

2  "  31/6r3(1/3)  " 


0.3^29 


where  r(l/3)  is  the  gamma  function  of  the  argument  l/3-  The  use  of 
the  relations  (10)  as  boundary  conditions'  along  the  sonic  line  reduces 
the  transonic  problem  of  the  flow  over  the  front  wedge  to  a  purely 
elliptic  problem  in- the  holograph  plane. 

In  addition  to  the  foregoing  conditions,  a  further  condition  is 
necessary  to  assure  that  the  solution  for  t  will  give  the  proper  scale 
when  transformed  to  the  physical  plane.  This  is  furnished,  for  example, 
by  the  following  expression  for  the  half -chord  of  the  profile,  found  by 
integrating  equation  (3a)  over  either  OB  or  0L: 
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2  f(y+i)(y0/a»-i) 

P*a*  L  2 


(Tl^±®w)  dTl 


If  the  chord  of  the  profile  is  given,  this  condition,  together  with  the 
previous  conditions  (7)  through  (10),  is  sufficient  to  determine  a  unique 
solution  to  the  problem. 

It  is  obvious  from  the  nature  of  the  boundary -value  problem  (and 
also  from  considerations  of  symmetry  in  the  physical  plane)  that  the 
solution  for  ijr  must  be  antisymmetric  with  respect  to  0.  The  problem 
can  be  simplified,  therefore,  by  discarding  the  lower  half  of  the  hodo- 
graph  and  replacing  it  by  the  condition 


t(T],0)  =  0  for  tj  1 


(12) 


The  resulting  problem  is  readily  solved  with  numerical  methods  by 
assuming ran  arbitrary  value  of  $  at  some  point  (as,  for  example, 
point  E),  solving  for  ijr  in  the  upper  half  of  the  hodograph  subject  to 
(the  conditions  (7),  (8),  (9) ,  (10),  and  (12),  and  then  adjusting  the 
(solution  to  satisfy  condition  (ll) . 


I  It  is  apparent  from  the  boundary  conditions  that  the  solution  of 
the  foregoing  problem  will  depend  on  the  value  of  the  parameter  0W, 

I  which  defines  the  position  of  the  upper  and  lower  boundaries  in  the 
I  hodograph.  This  parameter  is  related  to  the  more  familiar  transonic 
I  similarity  parameter  iQ  by  the  relation 

i  1’  . 


Mpg  -1  _  21/3 

[(7+1) (t/c )]2/3  " 


(13) 


where  t/c  is  the  thickness  ratio  of  the  complete  double-wedge  profile. 
In  references  2  and  3  "the  solution  of  the  foregoing  problem  has  been 
carried  out  for  four  values  of  9^, 


Ik 
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Problem  at  angle  of  attack.-  If  the  supersonic  region  at  the 
leading  edge  is  ignored,  the  boundary -value  problem  for  the  wedge  at 
angle  of  attack  appears  in  the  r ]f9  plane  as  shown  in  sketch  (d). 


6 


Sketch  (d) 

The  primary  difference  between  this  and  the  previous  sketch  is  that  the 
lines  OB  and  OL,  which  represent  the  surfaces  of  the  wedge,  have  each 
been  displaced  downward  by  an  amount  a,  where  a  is  the  angle  of 
attack  normalized  in  the  same  manner  as  the  other  angles  of  inclination 
(cf.  equation  (lb)).3 _ _  . .  ......  ,  ■ 

3In  reference  9>  Guderley  and  Yoshihara  find  it  convenient  to  obtain  the 
angle  of  attack  by  holding  the  profile  fixed  and  changing  the  inclina¬ 
tion  of  the  free  stream.  This  procedure,  if  applied  in  the  present 
case would  require  the  eventual  calculation  of  the  second  derivatives 
of  ijr  on  the  shock  polar.  The  present  procedure,  which  holds  the  free 
stream  fixed  and  changes  the  attitude  of  the  profile,  requires  the  cal¬ 
culation  of  only  a  first  derivative  of  at  the  surface  of  the  wedge. 
Since  the  accuracy  of  numerical  differentiation  decreases  with  increas¬ 
ing  order  of  the  derivative,  the  present  approach  is  to  be  preferred  in 
a  numerical  analysis. _ .  _ 
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The  -stream  function  at  angle  of  attack  will  be  denoted  here 
t>y  ♦=i|r(Tl»0;°O»  the  latter  notation  being  used  to  indicate  the  depend¬ 
ence  of  ^  upon  the  parameter  a.  The  function  must  satisfy  the 
differential  equation  (2),  which  is  now  written 

V)  -  2T]tee  =  0  (14) 

The  boundary  conditions  at  the  surface  of  the  wedge  now  require  that 

^(T],±0w-aja)  =  0  for  (15) 


while  the  condition  at  the  leading  edge  becomes 

^ - ->  0  for  tj - >-«,  -0w-cc <  0 <  0w-a 


(16) 


The  shock  polar  NAE  is  unaltered  from  the  previous  problem,  and  the 
condition  on  this  boundary  has  the  same  form  as  before.  The  conditions 
along  the  segments  BE  and  LN  of  the  sonic  line  are  now 


Vo,e;a)  + 


±Qy~Ct 


^e(°,e1;a) 

- -  d01  =  0 

[±(0!-e)f/3 


(17) 


where  the  lower  limit  of  the  integral  has  been  changed  in  accord  with 
the  displacement  of  the  points  B  and  L.  The  upper  signs  in  equation  (17) 
now  apply  for  1^0^0w-a  and  the  lower  signs  for  -0w-a^0-£-  1.  An 
expression  for  the  half -chord  of  the  profile  can  be  found  again  by  inte¬ 
grating  equation  (3a-)  over  the  line  OB  or  OL,  which  gives 


£ 

2 


2 

P*<% 


(7+l)(V0/a#-l) 

2 


o 

4'ife(i1#±0w-a;a)  dll 


(18) 


If  the  chord  of  the  airfoil  is  specified  -  say  the  same  as  at  zero  angle 
of  attack  -  then  the  foregoing  conditions  are  sufficient  to  determine  a 
solution.  No  simplification  based  on  symmetry  considerations  is  possible 
in  the  present  case. 

Perturbation  problem.  -  The  problem  of  the  preceding  section  conceiv¬ 
ably  could  be  solved  by  numerical  methods  -  though  with  great  labor  -  for 
arbitrary  values  of  a.  Efforts  in  this  direction  would  hardly  be  justi¬ 
fied,  however,  in  view  of  the  fundamental  omission  of  the  localized 
supersonic  flow  at  the  leading  edge.  It  is  more  reasonable  to  examine 
the  problem  for  vanishingly  small  a,  where  this  omission  is  valid  and 
where  there  is  hope  that  the  amount  of  labor  might  be  reduced. 
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To  proceed  along  these,  lines,  it  is  assumed  that  t(T),0ja)  may  he 
expanded  in  a  power  series  of  the  form 

y(ri,e;a)  =  *(-T),0;O)  +  a  ta(Tl,0>O)  +  0(a2) 


where,  for  present  purposes,  only  terms  to  order  a  need  he  retained. 

The  first  term  on  the  right  represents  the  solution  at  a  =  0  and  is 
thus  identical  with  the  function  Y(ti,0)  previously  introduced.  The 
second  term  will  he  abbreviated  hy  means  of  the  notation  '|r^>e)=1lra(Tl>0>0'' 
If  terms  of  0(a2)  are  discarded,  the  expression  for  ^  can  then 
he  written 

i|r(r),  0;a)  =  ^(ti, 0)  +  a  t*  (n,e) 

By  comparison  of  the  previous  boundary-value  problems  for  and  f,  a 
problem  for  the  perturbation  function  V  will  now  he  formulated. 

The  differential  equation  for  ’  follows  at  once  from  the  differ¬ 
ential  equations  (6)  and  (l4)  and  the  substitution  (19)-  I1:  is  obviously 

of  the  same  form  as  the  previous  equations,  that  is, 


^•2O'0r° 

The  boundary  conditions  appropriate  to  the  surface  of  the  wedge  are 
established  as  follows:  The  boundary  condition  (15)  for  +  is  first 
rewritten,  with  the  aid  of  the  substitution  (19) ?  in  form 


'F(Tb±0w-a)  +  a  (n>±0w-a)  =  0  (21) 

By  expanding  in  Taylor’s  series  about  the  lines  0  =  ±0W,  the  functions 
\j?  and  i|r’  can  be  written 

♦(q,±0w-a)  =  ^(ti,±0w)  -  a  $0(ti,±0w)  +  0(a2)  (22a) 

t,(Tb±ew"a)  =  ’Kb^w)  "  a  ^©(b^w)  +  °(a2)  (22b) 

If  these  expansions  are  substituted  into  equation  (2l)  and  +(11, ±0W)  set 

equal  to  zero  in  accord  with  the  boundary  condition  (7),  one  obtains 
finally  for  vanishingly  small  ,a 

V(t],±0w)  =  ^e(T],±0w)  for  rj^O  (23) 
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This  is  the  boundary  condition  for  V  appropriate  to  the  surface  of 
the  profile.  It  will  be  noted  that  the  condition  is  applied  in  the 
hodograph  at  the  original,  undeflected  location  of  the  surface 
(i.e. ,  6=±9V) .  The  condition  depends  for  its  application  on  a  knowledge 
of  the  basic  solution 

The  boundary  condition  for  at  the  leading  edge  follows 

directly  from  the  conditions  (8)  and  (l6).  It  is  the  same  as  the  corre¬ 
sponding  condition  for  if,  that  is, 

t* - for  T) — >-«,  -ev<6<6v  {2k) 

As  was  indicated,  the  functions  \|r  and  if  both  satisfy  the  same 
linear,  homogeneous  boundary  condition  on  the  shock  polar.  It  follows, 
as  in  the  case  of  the  differential  equation,  that  the  condition  for  ty* 
on  the  polar  is  again  the  same,  that  is, 

*'e  =  0  (25) 

for 

9  =  ±(l-T))'/l+T],  -I^ti^O 

.4 

The  treatment  of  the  boundary  condition  along  the  sonic  line  is 
complicated  by  the  fact  that  the  parameter  a  appears  in  the  condi¬ 
tion  (17)  as  a  term  in  the  lower  limit  of  the  integral.  For  simplicity, 
the  details  will  be  confined  here  to  the  upper  segment  BE  of  the  sonic 
line.  For  this  segment,  condition  (17)  becomes,  after  substitution 
from  equation  (19), 


*rj(0,e)  +  a ^(0,6)  +  k£f6  d0x  +  re(°M 


9w-a  (01  -0)2/3 


ew-a  (01 


2/3  d0l  = 0 


(26) 


applicable  for  1^0^0w-a.  To  simplify  this  equation,  the  first 
integral  is  rewritten 


/ 


0  ♦0(O,0l) 


^  -f 


*e(0,9i) 


ew-a  («1  -ef/3 


ae,  -/ 


0w-a  fyo,ex) 

T/ 3 


(01-0) 


(27) 


18 


MCA  IN  2832 


It  can  be  shown  from  Guderley’s  analysis  of  flow  at  a  convex  corner 
(reference  13)  that,  for  vanishingly  small  values  of  (0w-0),  the 
variation  of  if  along  the  sonic  line  must  be  of  the  form 

•• 

*(0,0)  =  c  (ew-e)4/3  (28) 

where  C  is  a  constant  for  any  given  value  of  0W.  Differentiating 
this  relation,  one  obtains 

1  /3 

tyo^)  ~  (0w-®i) 

Substitution  of  this  result  into  the  second  integral  of  equation  (27) 
yields  the  fact  that  this  integral  must  be  proportional  to  a4/3.  The 
first  integral  in  equation  (26)  can  thus  be  written 


,0  fe(0 ,ex) 


(01-0) 


275  d0l  = 


,0  *fl(O,0, ) 

e  ae,  +  0(a*/9) 

(8l-e)E/3 


The  second  integral  of  equation  (26)  can  be  treated  similarly  by  first 
rewriting  it  as 


re  re(  0,03.) 

V-a  (o^ef'3 


d0x  = 


(e^ef'3 


d0,  - 


To  deduce  the  variation  of  *’  for  vanishingly  small  (0w-0),  it  is 
first  noted  that  a  result  similar  to  equation  (28)  must  also  hold  for 
the  variation  of  *  relative  to  the  displaced  location  of  the  shoulder, 
that  is. 


*(0,0)  =  C  (0w-a-0)  ^ 


The  quantity  C  =  C(<x)  is  a  differentiable  function  of  ot  which  reduces 
to  C  when  a  =  0.  Since  a  will  eventually  be  made  less  than  any 
assignable  value  of  (0w-0),  this  expression  may  be  expanded  in  the  form 


*(.0,0)  =  C  (0w-0)' 


1  3  _8L_  .+  0(a2) 
.  3  ©w-0 
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Now  it  follows  from  the  definition  of  that 

ft  (O  0)  _  lim  ^(O>0)  -  *(0,0) 

CL  ^O  ~ 


Substitution  from  equations  (28)  and  (31)  thus  gives  for  the  variation 
of  f  in  the  vicinity  of  the  shoulder 


*»(o ,e) 


(ew-e)4/3  -h.c  (0w-e)1/3+  o(a) 


or 

♦*(0,0)  =  c»  (0w-©)4/3-  |  c  (ew-e)1/3 

where  C'=Ca(0).  This  means  that  for  vanishingly  small  (0w-0) 

♦*0  (0,®i)  -  (0w-0i)"2/3 

On  the  basis  of  this  result  equation  (30)  can  he  written 


(32) 


f 

0w-a 


**0(0,0!) 
(0x-0)2/3  ^ 


9  **0(0,61) 
(0^0)  2/3 


d0x  +  0(a1/3) 


(33) 


If  equations  (29)  and  (33)  are  substituted  into  equation  (26)  and  the 
boundary  condition  (10)  is  taken  into  account,  one  then  obtains  for 
vanishingly  small  a 


r6 


d0i  =  0 


(34) 


where  1^0^0W.  The  boundary  condition  for  f  along  the  upper  segment 
of  the  sonic  line  is  thus  the  same  as  the  condition  for  ijf.  The  same 
result  can  be  shown  to  hold  along  the  lower  segment. 

It  remains  to  impose  the  condition  that  the  chord  of  the  airfoil 
must  remain  unaltered  during  change  in  angle  of  attack.  To  express  this 
condition  in  terms  of  equation  (19)  is  first  substituted  into 

equation  (18)  to  obtain 


c  _  2 

'(7+1)  (vo/a*-l)~ 

V2  no 

r  -| 

2  1 

L  2  J 

J  ' 

-  OO 

(35) 


I 
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As  in  the  treatment  of  the  boundary  conditions  along  the  upper  and 
lower  boundaries,  Taylor's  expansion  gives 

1V'n,±Sw)-a  ^ee(ii,±©w)  +  o(a2)  (36a) 

+  o(a2)  (36b) 

It  can  be  inferred  directly  from  the  boundary-value  problem  for  f 
that  ^^(ri,  ±0W)  =  0,  so  that  the  term  involving  this  quantity  may  be 
dropped  from  equation  (36a).  Substitution  of  equations  (36)  into 
equation  (35)  and  application  of  the  previous  expression  (ll)  leads, 
for  vanishingly  small  a,  to  the  condition  that 

o 

T)  t'0(Tl,±ew)  dT}  =  0  (37) 

The  boundary  conditions  (23),  (24),  (25),  (34),  and  (37)  are  sufficient 
to  determine  the  solution  for  t '  in  the  hodograph. 

As  with  ijr,  the  boundary -value  problem  for  f  can  be  simplified 
from  considerations  of  symmetry.  Since  f  is  antisymmetric  with 
respect  to  0,  the  nonhomogene ous  boundary  condition  (23)  which  is 
imposed  on  f  along  the  upper  and  lower  boundaries  must  be  symmetric 
in  this  variable.  The  remaining  conditions,  which  are  all  homogeneous, 
are  also  symmetric.  It  follows  that  itself  must  be  a  symmetric' 

function  of  0.4  The  problem  can  therefore  be  simplified  by  again 
eliminating  the  lower  half  of  the  hodograph  and  substituting  in  this 
case  the  condition  that 

t'^TbO)  =  0  for  Tl^“  1  (38) 


4 This  result  can  also  be  argued  directly  from  considerations  in  the 
physical  plane.  It  is  necessary  to  make  two  observations  as  follows: 
(l)  Since  the  profile  itself  is  symmetric  about  the  chord  line,  the 
flow  field  at  a  negative  angle  of  attack  must  be  the  inverted  image 
of  the  flow  field  at  an  equal  positive  angle.  (2)  To  be  consistent 
with  the  basic  perturbation  assumption,  it  must  be  presumed  that  all 
changes  in  the  flow  field  are  smooth  functions  of  angle  of  attack 
at  a  =  0.  These  statements  taken  together  imply  that  the  vertical 
distance  between  any  two  points  of  equal  tj  and  corresponding  posi¬ 
tive  and  negative  0  is,  to  a  first  order,  unaffected  by  angle  of 
attack.  It  follows  that,  for  sufficiently. small  a,  the  increments 
in  ^  and  t  between  the  two  points  are  equal  and  hence,  on  the  basis 
of  equation  (19) ,  that  the  value  of  V  at  the  two  points  is  the  same. 


^0-(n»*0w-a)  = 

*e(Tb±0w-a)  = 
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The  problem  which  is  finally  to  be  solved  is  thus  as  summarized  in 
sketch  (e).  The  boundaries  for  this  problem  are  identical  with  those 
used  to  obtain  if.  The  boundary  conditions  are  also  identical  insofar 
as  the  shock  polar  and  sonic  line  are  concerned.  The  only  differences 
between  the  two  problems  are  in  the  conditions  imposed  along  the  bound¬ 
aries  OB  and  OA.  As  was  the  case  with  if,  the  solution  here  must  be  a 
function  of  9V. 


0 


Because  of  the  nature  of  the  integral  condition  along  the  upper 
boundary  OB,  a  direct  solution  for  _♦*  is  not  feasible  by  numerical 
methods.  To  obtain  a  solution,  therefore,  the  problem  is  broken  down 
into  two  subsidiary  problems  by  means  of  the  substitution 


'  f  =  *Y+  b*’B  (39) 

where  b  is  a  constant  whose  value  is  to  be  determined.  Boundary -value 
problems  for  and  ijf'g  are  then  defined  as  shown  in  sketch  (f ) . 
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In  both  problems  the  integral  condition  along  the  upper  boundary  OB  is 
ignored  for  the  time  being,  and  an  arbitrary  condition  is  introduced 
instead  at  the  point  E.  In  the  problem  for  t’A,  only  the  nonhomoge- 
neous  condition  (23)  is  imposed  at  the  upper  boundary,  and  the  condition 
at  E  is  the  homogeneous  one  that  =0.  In  the  problem  for 

the  homogeneous  condition  =0  is  imposed  along  OB,  and  the  condi¬ 
tion  at  E  is  that  has  an  arbitrary  value  ^  0.  The  conditions 

at  the  remaining  boundaries  are  the  same  as  in  sketch  (e)  and  are  there¬ 
fore  not  repeated  here.  It  is  apparent  that  a  superposition  of  ty'A 
and  will  constitute  a  solution  of  the  original  problem  provided  the 

value  of  b  is  adjusted  so  that  the  integral  condition  (37)  is  satisfied 
on  the  upper  boundary.  The  necessary  equation  for  b  is  found  by  sub¬ 
stituting  the  expression  (39)  into  condition  (37)  and  is 


(t],0w)  dr) 


/ 


o 

■n  Vbq^j^w) 


(bo) 


Relations  for  quantities  in  physical  plane.-  To  complete  the  funda¬ 
mental  analysis,  relations  must  be  established  between  and  the 

relevant  quantities  in  the  physical  plane.  Let  x  =  x(t|,0)  and  y  =  y(T],0) 
denote  the  coordinates  at  which  a  given  velocity  T],0  is  found  in  the 
physical  plane  when  the  profile  is  at  zero  angle  of  attack.  As  shown  in 
reference  3  (pp.  29-31),  the  transformation  equations  ( 3) ,  when  applied 
to  the  case  of  zero  angle  of  attack  (and  written  in  the  present  notation), 
can  be  put  in  the  dimensionless  form 


(2T)f0dT  +  t^d©) 


(kla) 


[(7+i)(t/c)]1/3  d(j)  =  (2^-/3 ( +  *ede)  = 


(4lb) 


where  Iw  represents  the  integral 


Tltfl('n70w)  dTl 
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By  taking  the  origin  of  the  physical  coordinates  at  the  leading  edge  and 
introducing  the  notation  X  =  x/c  and  Y  =  [ (7+1) (t/c)]1/3(y/c)  equa¬ 
tions  (4l)  can  he  integrated  to  give 

'^,6)  =  ~  J  (2Tj^0dT)  +  ^de)  (43a) 

W  C 

(pa 

Y(n,0)  =  v  (43b) 

4lw 

The  integration  in  equation  (43a)  is  performed  in  the  hodograph  over  any 
curve  C  which  begins  at  3  =  -00  and  ends  at  the  point  4,0.  The 
generalized  coordinates  X  and  Y  at  which  the  same  velocity  r\>6  is  to 
be  found  when  the  airfoil  is  at  angle  of  attack  are  given  correspondingly 
■by 


^(4*®#®)  =  -j -  J  (2T)\(f0dT)  +  \|r  d0) 

.  4IW  Jq  ^ 

(44a) 

*».,«>-  (29»)l/\ 

4iw 

(44b) 

The  integration  in  equation  (44a)  is  considered  to  be  taken  over  the 
same  curve  C  as  before.5  The  integral  Iw  is  now  given  by 

Iw  —  h  ( 4  >  ^w-01)  3-4 

m 

'  —  uu 

It  can  be  shown  from  equations  (19),  (36),  and  (37)  that  for  vanishingly 
small  a  . 

.H* 

II 

(^6) 

Equations  (44)  can  now  be  specialized'  in  the  light  of  the 
perturbation  assumption.  This  assumption  implies  at,  once  that 
dinates  X  and  Y  in  the  physical  plane  must  be  expressible  in 

basic 
the  coor- 
the  form 

If.  c  Hes  slightly  outside  the  domain  in  which  f  is  defined  -  as 
will  be  the  ease,  for  example,  when  the  integration  is  taken  over  the 
upper  surface  of  the  wedge  in  its  undisplaced  position  -  is  to  be 
thought  of  as  being  continued  analytically  outside  the  boundary. 
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X(ti,0;cc)  =  X(q,0)  +  a  X*  (tj,©)  (47a) 

f  . 

Y(n,0;a)  =  Y(ri,e)  +  a  Y»  (T!,0)  (47b) 

where  X’ (t),0)  HXa(T],0;O)  and  Y*  (t],0)  =  Ya(rj.0;O) .  If  expressions  (19) 
and  (47)  are  substituted  into  equations  (44),  and  equations  (43)  and  (46) 
are  taken  into  account,  the  following  relations  are  finally  obtained 
for  X*  and  Y*  in  terms  of  \|r' : 


(48a) 


(48b) 


The  foregoing  equations  (48)  give  the  initial  rate  of  movement  with 
angle  of  attack  of  a  point  of  fixed  velocity  r|,0.  One  requires  for 
practical  application,  however,  the  rate  of  change  of  t]  and  0  at  a 
point  of  fixed  location  X,  Y.  Equations  relating  the  two  sets  of  deriv¬ 
atives  can  be  obtained  as  follows:  If  tj  and  0  are  regarded  in  the 
physical  plane  as  functions  of  X,  Y,  and, a  -  that  is,  tj  =  T](X,  Y;  a) 
and  0  =  0(X,  Y;  a)  -  then  the  corresponding  total  differentials  are 

dq  =  T)xdX  +  TjydY  +  T|ada,  (49a) 

d0  =  0xdX  +  0ydY  +  0ada  (49b) 

Consistent  with  the  basic  perturbation  assumption,  r)  and  0  can  be 
written 

T](XjY,a)  =  n(X,Y)  +  a  t!»(X,y)  (50a) 

0(X,Y,a)  =  0(X,Y)  +  a  0’ (X,Y)  (50b) 

where  t}  and  0  represent  the  conditions  at  a  given  point  X,  Y  at  zero 
angle  of  attack  and  tj*  and  0'  are  defined  by  ri,(X,  Y)  =  ria(X,  Y;  0) , 
0’(X,  Y)=0a(X,  Yj  0).  In  view  of  equations  ( 50 ),  equations  (49)  can 
be  written  for  vanishingly  small  a 

dt]  =  fj^cLX  +  fJydY  +  T]rda  -  (51a) 

■  .  ■  ;  - 

d0  =  e^dX  +  0ydY  +  0'da  (51h) 


X*(ti,0) 

Y»(ti,0) 


ji-  f  (2^’ed,  +  ye) 
(aur 


,1/3 


4i 


r 


w 


4U 
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Similarly,  from  equations  (47),  one  can  write  for  the  differentials 
of  X  and  Y  as  functions  of  t],0,  and  a 


dX  =  X^dT)  +  XgdS  +  X*da 
•  dY  =  Y^dT]  +  Ygd0  +  Y*da 

from  which 

X^d-q  +  Xgd0  =  dX  -  X'da  (52a) 

Y^dr]  +  Yed0  .  dY  -  Y'da  (52b) 

Solution  of  equations  (52)  for  dr)  and  d0  and  comparison  of  the  results 
with  the  alternative  expressions  (5l)  gives  finally  for  t)*  and  0' 


V'  -  XeY' 

-  XgY, 


(53a) 


V*  -  X^Y* 

X^Yg  -  XgYr, 


(53b) 


These  equations  can  he  put  in  more  directly  useful  form  by  evaluating 
the  derivatives  of  X  and  Y  from  equations  (43)  and  substituting  for 
X*  and  Y'  from  equations  (48).  There  results  finally 


V(x,Y) 


(2r}t  gdt]  +  i^'^d©) 


(54a) 


e*(x,Y) 


+t’ d0)  -2^+' 


(54b) 


By  means  of  these  equations  the  initial  rate  of  change  of  t)  and  0  at 
some  fixed  point  in  the  physical  plane  can  be  calculated  corresponding 
to  any  chosen  location  in  the  holograph.  The  coordinates  at  which  these 
derivatives  apply  are  found  from  the  solution  at  zero  angle  of  attack 
by  means  of  equations  (43) . 


2  6 
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The  foregoing  equations  are  considerably  simplified  when  applied  at 
the  surface  of  a  wedge  profile.  Here  the  boundary  condition  is  that^ 
is  constant  on  a  line  of  constant  6  (cf.  equation  (7) ) ,  with  the  result 
that  4^=0.  Equation  (54a),  for  example,  can  thus  be  written  as 
simply 

■n 

V(x,±o)  - - - -  f  w)  dTi  (55) 

■ii+e(Tl,±0w) 


where  the  upper  signs  pertain  to  the  upper  surface  and  the  lower  signs 
to  the  lower  surface.  The  corresponding  rate  of  change  of  pressure 
coefficient  is  found  by  differentiating  equation  (4)  with  respect  to 
angle  of  attack.  If  a  is  used  now  to  denote  the  true  angle  of  attack 
(related  to  the  previously  used,  normalized  angle  of  attack  by  an  equa¬ 
tion  like  (lb)),  such  differentiation  then  gives 


1 _ 

(To  /a»-l)1/2 


Here  tj'  is  still  the  derivative  with  respect  to  the  normalized  angle 
as  given  by  equations  (54a)  or  (55)-  With  the  aid  of  equations  (5)  and 
(13),  this  result  can  be  rewritten 

[(7+l)(t/c)]1/3  =  "  2(2@w)  /  V  (56) 

x  acL  xa=o 

It  can  be  seen  from  equation  (55)  and  the  symmetry  properties  of  if 
and  ijr1  that  T)’  must  be  of  equal  magnitude  but  opposite  sign  on  the 
upper  and  lower  surfaces  of  the  profile.  If  the  local  lift  coefficient 
is  represented  by  Ap/lo  -  (Plower  "  PUpper)/(lo>  it;  then  Allows  from 
equation  (56)  that 


[(r+iXt/c)]1'3 


da 


=  4(20W)1/3V  (X,+0) 


-■a=o 


(57) 


where  the  notation  T|*(X,+0)  indicates  that  the  value  is  to  be  taken  on 
the  upper  surface  of  the  profile.  Substitution  from  equation  (55)  gives 
finally 


1(7+1)  (t/c)]V3 


jjAp/qp) 

da 


4(2@w)1/3 


Ja=o 


n+0(^i,0w)  . 


P 

J  ■nt^T),©^)  dq  (58) 


By  means  of  this  equation,  the  initial  rate_  of  growth  of  lift  at  any 
chordwise  station  can  be  obtained.  Since  if  and  are  both  functions 
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of  the  parameter  0W,  the  generalized  quantity  which  appears  on  the  left- 
hand  side  of  equation  (57)  is  also  a  function  of  this  parameter.  These 
results  are  in  conformity  with  the  rules  for  transonic  similarity  (see, 
for  example,  reference  11) . 


METHOD  OF  SOLUTION 


As  in  the  previous  calculations  of  ijjr,  the  "boundary-value  problems 
for  and  i|f’g  can  he  solved  through  the  use  of  finite -difference 

equations  and  relaxation  techniques.  A  detailed  description  of  the 
general  method  has  been  given  in  reference  3  and  need  not  be  repeated 
here.  Most  of  the  necessary  finite -difference  equations  -  notably  the 
tedious  ones  along  the  shock  polar  and  sonic  line  -  can  be  taken  over 
directly  from  the  previous  work.  The  only  equations  which  need  be 
altered  are  those  directly  influenced  by  the  change  in  boundary  condi¬ 
tions  on  the  upper  boundary  and  on  the  horizontal  axis.  The  only  real 
difficulty  from  this  source  is  encountered  in  the  solution  for  \Jr*^  in 
the  vicinity  of  the  shoulder  (point  B  in  sketch  (f)).  At  the  shoulder 
itself,  the  boundary  conditions  require  a  singularity  in  the  first 
derivatives  of  +  which  means  that  any  purely  numerical  treatment 
would  be  of  doubtful  validity  in  the  vicinity  of  this  point.  This  dif¬ 
ficulty  is  overcome  by  subtracting  out  an  analytical  solution  of  the 
proper  singular  form  and  then  working  locally  with  the  difference 
between  this  solution  and  the  desired  unknown.  The  singular  solution 
is  obtained  from  the  general  results  of  Guderley  (reference  13)  and  is 
expressed  in  terms  of  hypergeometric  functions.  The  details  of  this 
and  other  matters  regarding  the  numerical  calculations  for  the  front 
half  of  the  profile  are  given  in  appendix  A. 

With  the  solution  known  for  the  front  half  of  the  profile,  the 
calculation  of  the  lift  on  the  rear  half  is  a  simple  matter.  The  com¬ 
putations  are  carried  out  in  the  physical  plane  and  are  based  on  the 
characteristics  net  previously  constructed  for  the  flow  over  the  rear 
wedge  at  zero  angle  of  attack  (see,  for  example,  fig.  k  of  reference  3). 
Starting  from  the  known  solution  for  "'If*,  one  first  employs  equations  (48) 
to  compute  the  initial  rate \ of  movement  of  the  points  at  which  the  Mach 
lines  of  the  basic  characteristics  net  meet  the  sonic  line.  Using  these 
results  and  the  known  slope  of  the  segments  of  the  basic  net,  one  then 
proceeds  stepwise  along  consecutive  downgoing  Mach  lines,  calculating 
the  initial  rate  of  movement  of  successive  intersection  points  on  each 
line.  By  application  of  the  proper  boundary  conditions  at  the  surface 
of  the  wedge,  the  value  of  tj*  at  the  surface  is  finally  determined,  and 
from  this  the  initial  distribution  of  lift  is  calculated.  The  details 
of  the  procedure  are  given  in  appendix  B. 
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RESULTS  AED  DISCUSSION 


Calculations  of  the  lift  have  been  carried  out,  following  the 
methods  just  outlined,  for  the  same  values  of  0W  used  in  the  work 
at  zero  lift,  namely,  1.3,  1.6,  2.4 ,  and  4.2.  These  values  correspond, 
respectively  (see  equation  (13) ),  to  values  of  the  similarity  param¬ 
eter  £0  of  1.058,  0.921,  0.703,  and  0.484. 

To  illustrate  the  results  for  the  front  wedge  in  the  hodograph, 
figures  1  to  3  have  been  prepared  showing  the  variation  of 

and  x|/ 1  for  0W  =  I.6.6  The  results  for  iJM-g  (fig.  l)  are  only  slightly 

different  from  those  previously  shown  for  if  in  figure  3  of  reference  3* 
As  before,  a  rapid  (but  regular)  variation  is  apparent  in  the  dependent 
variable  in  the  vicinity  of  the  point  q  =  0,  0  =  1.  The  results  for 
(fig.  2)  show  a  rapid  variation  near  the  point  q  =  0,  0  =  0W.  This  is 
a  consequence  of  the  previously  mentioned  singularity  in  the  first 
derivatives  of  at  that  point.  The  values  of  if’  (fig.  3)  are 

found  in  the  present  case’ from  the  equation  hf*  =  -  0-5348  if'g 

(cf.  equation  (39))*  They  exhibit  the  same  behavior  as  does  in 

the  vicinity  of  the  singular  point  but  differ  markedly  in  other  parts  of 
the  field.  For  reference,  the  numerical  values  from  which  figures  1  and  2 
were  plotted  are  given  in  tabular  form  at  the  end  of  the  report. 

The  complete  results  for  the  lift  of  the  profile  are  given  in 

figures  4  through  7*  These  results  will  be  discussed  in  the  following 
paragraphs. 


Chordwise  Distribution  of  Lift 


Figure  4  is  a  plot  of  the  calculated  lift  distribution,  in  tran¬ 
sonic  similarity  form,  for  the  four  values  of  considered  in  the 

present  work.  Also  shown  are  the  results  for  |Q  =  0  (Mq  =  l)  given  by 
Guderley  and  Yoshihara  in  reference  9*  It  is  convenient  for  purposes 
of  discussion  to  think  of  a  similarity  plot,  such  as  that  of  figure  4, 
as  pertaining  to  fixed  values  of  t/c  and  7.  From  this  standpoint,  an 
increase  from  zero  in  the  similarity  parameter  |D  is  equivalent  to  an 
increase  from  1  in  the  free -stream  Mach  number  Mq.  For  simplicity, 
this  point  of  view  will  sometimes  be  adopted  in  the  descriptions  that 
follow. 


6For  the  calculation  of  ty’-g  in  this  example,  use  was  made  of  236 

lattice  points  distributed  as  shown  for  f  in  figure  2  of  reference  3- 
For  380  points  were  used  with  a  distribution  appropriate  to  the 

altered  behavior  of  the  dependent  variable. 
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The  lift  distributions  of  figure  4  are  all  of  the  same  general 
shape.  In  all  cases  the  calculated  lift  tends  toward  infinity  at  the 
leading  edge  of  the  profile.  This  type  of  Iresult,  which  is  of  course 
physically  impossible,  is  well  known  from  the  linear  theory  of  airfoils 
at  subsonic  speeds.  It  is  a  result  of  the  obvious  failure  of  the  small - 
disturbance  approximations  to  conform  with  the  actual  phenomena  in  the 
vicinity  of  the  leading  edge.  This  local  failure  of  the  theory  is 
known  in  the  linear,  subsonic  case  to  be  of  little  consequence  insofar 
as  the  over-all  lift  is  concerned.  It  may  be  presumed  that  a  similar 
situation  exists  here. 

As  one  proceeds  rearward  from  the  leading  edge,  the  lift  distribu¬ 
tion  falls  more  or  less  rapidly,  reaching  a  value  of  zero  directly 
forward  of  the  shoulder.  This  latter  result  could  have  been  foreseen, 
since  the  speed  on  both  the  upper  and  lower  surfaces  has  a  fixed 
(i.e.,  sonic)  value  at  this  location.  Directly  to  the  rear  of  the 
shoulder,  the  lift  distribution  starts  anew  from  zero.  This  must 
obviously  be  the  case,  since  the  expansion  from  sonic  speed  is,  in 
Prandtl -Meyer  flow,  a  unique  function  of  the  local  turning  angle,  which 
is  the  same  for  both  surfaces.  Rearward  from  the  shoulder  the  lift 
increases  monotonically  to  a  relatively  small,  finite  value  at  the 
trailing  edge. 

Over  the  front  wedge,  the  four  curves  of  the  present  study  exhibit 
a  uniform  progression  with  respect  to  £Q.  The  curve  of  Guderley  and 
Yoshihara,  however,  crosses  the  present  curves  at  several  points.  The 
reasons  for  this  are  not  clear,  though  it  is  highly  unlikely  that  such 
a  result  could  be  in  fact  correct.  The  observed  behavior  may  be  due  to 
some  consistent  inaccuracy  in  the  present  numerical  approach  or  to  the 
approximations  introduced  by  Guderley  and  Yoshihara  in  satisfying  the 
boundary  conditions  for  the  interdependent  portion  of  the  supersonic 
expansion  fan.  Over  the  rear  wedge,  the  present  computations  give  vir¬ 
tually  a  single  curve  for  the  four  values  of  |Q.  There  is  again,  how¬ 
ever,  a  small  inconsistency  with  the  results  given  by  Guderley  and 
Yoshihara.  This  is  as  might  be  expected  if  the  calculated  flow  over  the 
front  wedge  is  in  error  in  either  case. 


Lift -Curve  Slope 


Figure  5  shows  the  generalized  slope  of  the  lift  curve  at  zero 
angle  of  attack  plotted  as  a  function  of  the  transonic  similarity  param¬ 
eter.  Results  obtained  on  the  basis  of  the  transonic  small -disturbance 
theory  are  shown  by  three  solid-line  curves.  Each  of  these  curves  con¬ 
sists  of  two  segments  separated  by  a  gap  within  which  the  curve  cannot 
be  defined  on  the  basis  of  the  available  results.  The  uppermost  of  the 
three  curves  gives  the  lift  of  the  complete  profile;  the  other  two  show 
the  division  of  lift  between  the  front  and  rear  wedges. 


30 


NACA  TN  2832 


The  left-hand  segment  of  each  of  the  curves  in  figure  5  shows  the 
variation  of  lift-curve  slope  over  most  of  the  range  of  flight  speed  in 
which  the  bow  wave  is  detached,  which  is  the  range  of  primary  concern 
in  the  present  analysis.  The  calculated  points  from  which  these  curves 
were  drawn  are  shown  in  the  figure.  The  points  denoted  by  squares  were 
obtained  by  mechanical  integration  of  the  lift-distribution  curves  of 
figure  k. 7  The  circled  points  on  the  vertical  axis  were  located  on  the 
basis  of  the  work  of  G-uderley  and  Yoshihara. 

The  right-hand  segment  of  the'  curves  in  figure  5  shows  the  varia¬ 
tion  of  lift-curve  slope  in  the  range  of  flight  speed  in  which  the  bow 
wave  is  attached  and  the  flow  is  completely  supersonic.  To  the  order 
of  accuracy  of  the  present  theory,  this  condition  exists  for  the  double¬ 
wedge  profile  at  zero  angle  of  attack  when  £0  ^  21/^3  =  1.2 60.8  Above 
this  value,  results  completely  consistent  with  the  fundamental  assump¬ 
tions  of  the  transonic  small-disturbance  theory  can  easily  be  obtained 
by  analytical  methods.  To  this  end,  one  need  only  presume  that  the 
speed  is  constant  on  each  straight-line  portion  of  the  airfoil  surface, 
a  condition  which  is  actually  fulfilled  over  most  of  the  pertinent 
range  of  ^Q.  The  necessary  procedures  are  outlined  in  appendix  C.  To 
the  accuracy  of  the  transonic  small-disturbance  theory,  the  results 
provide  an  exact  solution  for  the  lift-curve  slope  of  the  front  wedge 
for  all  values  of  £Q  in  the  range  of  completely  supersonic  flow.  For 
the  rear  wedge  -  and  hence  for  the  complete  profile  -  the  solution  is 
exact  down  to  a  limiting  value  of  fjQ  somewhat  greater  than  1.260. 

Below  this  limit  the  interaction  of  the  shock  wave  from  the  bow  and  the 
expansion  fan  from  the  shoulder  influences  the  flow  over  the  rear  wedge, 
with  the  result  that  the  condition  of  constant  speed  is  not  satisfied. 
The  position  of  this  limit  is  difficult  to  determine  exactly.  As  shown 
in  appendix  C,  however,  it  must  lie  at  a  value  of  less  than  I.287. 

The  curves  for  the  rear  wedge  and  complete  profile  are  thus  approximate 
for  at  least  a  portion  of  the  interval  from  1.287  to  1.260  and  are 
therefore  shown  dotted  in  this  range.  It  can  be  demonstrated  that 
inclusion  of  the  interaction  effects  in  the  analysis  would  cause  an 
increase  in  the  computed  lift  for  the  rear  wedge.  Exact  results  would 
thus  lie  somewhere  above  the  dotted  portion  of  the  curves  in  figure  5- 


7As  in  the  earlier  calculations  of  the  drag  coefficient  at  zero  angle 
(cf.  page  36  of  reference  3),  the  integration  over  a  small  interval 
near  the  leading  edge  was  carried  out  analytically  on  the  basis  of  an 
asymptotic  representation  of  the  solution  in  the  holograph  plane. 

8 Attachment  of  the  wave  takes  place  at  the  somewhat  lower  value  of 

i0  =  3/(b)2/3  =  1.191.  For  1.191  <  <1.260  the  wave  is  attached 

but  the  flow  behind  it  is  still  subsonic. 
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The  most  interesting  aspect  of  figure  5  is  the  behavior  of  the  lift 
in  the  vicinity  of  shock  attachment.  Despite  the  gap  in  the  curves  in 
this  vicinity,  i't  is  obvious  that  the  lift-curve  slope  of  the  complete 
profile  must  attain  a  maximum  somewhere  in  the  range  from  =  I.O58 
to  |0  =  I.287.  This  is  in  marked  -  and  somewhat  surprising  -  contrast 
to  the  previous  results  for  the  drag  coefficient  at  zero  angle  of  attack, 
which  was  found  (reference  2)  to  decrease  monotonically  as  the  similarity 
parameter  increased  above  zero.  The  peak  in  the  curve  in  the  present 
case  is  accompanied  by  a  similar  variation  in  the  lift-curve  slope  of 
the  front  wedge.  The  results  for  the  rear  wedge  may  or  may  not  pass 
through  a  minimum  in  the  same  range  of  £0. 

A  determination  of  the  exact  shape  of  the  curves  in  the  vicinity  of 
shock  attachment  is  not  feasible  on  the  basis  of  the  present  laborious 
methods.  The  existing. curve  for  the  complete  profile  does  show  a  maximum 
in  the  range  of  completely  supersonic  flow,  but  this  is  in  the  portion 
of  the  range  in  which  the  computed  curve  is  known  to  be  erroneously  low. 
If  exact  results  were  available  for  all  values  of  £  ,  the  maximum  would 
undoubtedly  be  somewhat  higher  and  displaced  somewhat  to  the  left.  The 
infinity  which  appears  in  the  slope  of  the  curve  at  £Q  =  1.260  (see 
appendix  C)  would  probably  disappear  as  well.  The  lift  of  the  rear 
wedge,  which  now  goes  to  zero  at  £0  =  1.260,  would  presumably  remain 
finite  throughout. 

Within  the  transonic  range  itself,  the  curves  of  figure  5  show 
little  variation  for  some  distance  above  a  similarity  parameter  of  zero. 
This  is  in  accord  with  Guderley's  recent  analytical  study  of  two- 
dimensional.  flows  with  a  free-stream  Mach  number  close  to  1  (refer¬ 
ence  l4) .  The  results  of  Guderley’s  work  imply  that,  to  the  acchracy 
of  the  small -disturbance  theory,  the  lift-curve  slope  does  not  vary  as 
the  free-stream  Mach  number  passes  .through  1.  The  curves  of  figure  5 
have  been  faired  so  as  to  conform  with  this  requirement.  It  is  apparent 
that  Guderley's  result  of  zero  variation  may  be  taken  as  a  good  working 
approximation  even  at  Mach  numbers  some  distance  removed  from  1.  The 
same  result  was  found  in  reference  3  vith  regard  to  the  drag  coefficient 
of  the  complete  profile  at  zero  angle  of  attack. 

Over  most  of  the  range  of  completely  supersonic  flow,  the  lift- 
curve  slope  of  the  complete  profile  exhibits  the  type  of  variation  well 
known  from  linear  theory.  This  latter  theory  gives  for  the  lift-curve 
slope  of  all  thin  profiles 


Ac 1  .  4  , 


(59) 


\ 


32 


MCA  TN  2832 


which  can  be  written  in  terms  of  the  transonic  similarity  variables  as 


1(7+1)  (t/c)]1/3  ^  =  ~ T/T  (60) 

^O 


The  dashed  curve  in  figure  5  is  based  on  this  equation.  There  is  con¬ 
siderable  quantitative  difference  between  the  linear  and  nonlinear 
results  for  values  of  £Q  just -above  I.287.  As  gQ  increases,  how¬ 
ever,  the  curves  given  by  the  two  theories  appear  to  converge.  This 
latter  behavior  is  in  accord  with  'Spre  iter's  considerations  (refer¬ 
ence  11)  regarding  the  basic  relationship  between  the  linear  and  non¬ 
linear  theories. 

To  put  the  results  in  more  familiar  form,  the  lift-curve  slope  of 
the  complete  profile  has  been  replotted  in  figure  6  as  a  function  of 
Mach  number  for  7  =  1.4.  The  results  of  linear  theory  give  a  unique 
curve  defined  by  equation  (59)-  The  nonlinear,  transonic  theory  pro¬ 
vides  a  family  of  curves  with  thickness  ratio  as  a  parameter.9  As 
would  be  expected,  the  range  of  Mach  numbers  over  which  the  linear 
theory  is  a  poor  approximation  becomes  smaller  as  the  thickness  ratio 
is  reduced.  It  can  be  reasoned,  in  fact^  that  the  nonlinear  results 
must  tend  toward  the  results  of  the  linear  theory  as  t/c — >  0. 


Center  of  Lift 


Figure  7  shows  the  chordwise  position  of  the  center  of  lift  (x/c)z 
as  a  function  of  the  transonic  similarity  parameter.  The  arrangement 
of ^ the  figure  parallels  that  of  figure  5.  As  before,  the  indicated 
points  were  calculated  on  the  basis  of  the  lift  distributions  of 
figure  4.  The  curve  in  the  range  of  completely  supersonic  flow 
(  £0  =l*26o)  was  obtained  by  means  of-  the  equations  of  appendix  C. 

Only  results  for  the  complete  profile  are  shown...  .  ..  .. 


Here,  as  in  reference  2,  a  multiplicity  of  curves  could  be  obtained 
for  each  thickness  ratio  by  using  expressions  for  the  pressure  coef¬ 
ficient  and  similarity  parameter  different  from  those  of  equations  (4) 
and  (13) .  In  view  of  the  recent  developments  outlined  in  footnote  2, 
such  complications  now  appear  to  be  of  lessened  significance. 


MCA  IN  2832 


33 


The  movement  of  the  center  of  lift  with  increasing  Mach  number  is 
of  some  interest.  At  a  free-stream  Mach  number  of  1,  the  results  of 
Guderley  and  Yoshihara  indicate  a  position  about  29  percent  of  the 
chord  aft  of  the . leading  edge .  As  the  Mach  number  is  increased;  the 
center  of  lift  first  moves  forward,  slowly  in  the  initial  stages  and 
then  more  rapidly  as  the  condition  for  shock  attachment  is  approached. 
In  the  completely  supersonic  range,  this  trend  is  reversed;  the  center 
of  lift  then  moves  aft  toward  the  midchord  location  given  by  linear 
theory.  Apparently,  the  reversal  of  the  direction  of  motion  must  take 
place  rather  suddenly  in  the  vicinity  of  shock  attachment.  The  limit 
of  forward  movement  cannot  be  specified,  except  to  say  that  it  must  lie 
somewhere  ahead  of  22  percent  of  the  chord  (and  probably  aft  of  the 
leading  edge).  The  dotted  (i.e.,  inexact  )  portion  of  the  curve  passes 
precisely  through  the  quarter-chord  point  at  £Q  =  1.260.  (The  corre¬ 
sponding  lift  distribution  is  one  of  uniform  lift  on  the  front  wedge 
and  zero  lift  on  the  rear.)  Because  of  the  interaction  effects 

previously  discussed,  an  exact  result  would  lie  somewhat  above  the 
dotted  curve .  > 


CONCLUDING  REMARKS 


The  present  calculations  add  support  to  the  growing  conclusion 
(see  references  2,  6,  7,  8,  and  14)  that  no  marked  changes  take  place 
in  characteristics  of  airfoil  sections  as  the  free-stream  Mach  number 
passes  through  1.  The  establishment  of  this  conclusion  must  be  regarded, 
in  fact,  as  one  of  the  major  successes  of  recent  research  in  transonic 
flow.  In  the  present  case,  as  in  the  previous  study  of  the  drag  coef- 
lcient  at  zero  lift,  the  variation  of  the  aerodynamic  quantities  with 
free-stream  Mach  number  is  most  rapid  in  the  vicinity  of  shock  attach- 
ment.  Unlike  the  behavior  of  the  drag  coefficient,  however,  the  varia- 
ions  here  are  large  and  characterized  by  a  sudden  reversal  in  the  sign 
01  the  derivative.  In  drawing  conclusions  from  these  results  it  must 
be  remembered,  of  course,  that  the  theory  assumes  an  inviscid  medium 
and  an  airfoil  of  small  thickness  and  infinite  span.  It  also  assumes, 
in  effect,  that  the  angle  of  attack  Is  of  an  order  smaller  than  the 
thickness  ratio.  To  what  extent  the  results  will  be  valid  for  viscous 
flows  about  finite -span  airfoils  at  practically  usable  values  of  the 
thickness  ratio  and  angle  of  attack  is  difficult  to  say.  The  effects  of 
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finite  span,  for  example,  will  surely  cause  a  reduction  in  the  variations 
near  shock  attachment.  In  the  present  state  of  theoretical  development, 
the  study  of  these  effects  is  a  task  for  experiment. 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  Aug.  1,  1952 
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APPENDIX  A 


SOLUTION  OF  BOUNDARY-VALUE  PROBLEM 
FOR  FRONT  WEDGE  IN  HODOGRAPH  PLANE 


'  The  solution  of  the  boundary-value  problems  for  and  V'g  was 

accomplished  by  finite-difference  methods  similar  to  those  developed  for 
the  calculation  of  V  in  reference  3»  The  description  here  will  be 
limited  to  the  few  features  wherein  the  present  work  departs  from  that 
discussed  in  the  earlier  paper.  (See  general  remarks  under  METHOD  OF 
SOLUTION . )  The  notation  and  sketches  follow  the  conventions  used  in  * 
reference  3* 


Finite -Difference  Equations  Common  to  Both  Problems 


The  only  finite-difference  equations  common  to  the  problems  for 
an^  ^'b  but  n°t  found  in  the  problem  for  V  derive  from  the  bound¬ 
ary  condition  on  the  horizontal  axis  (see  sketches  (e)  and  (f)).  This 
condition  is  given  for  both  problems  by  equation  (38)  and  is  ^  'g  (tj,0)  =  0 
for  T)  5--1.  In  the  previous  work,  the  finite-difference  equations  for 
lattice  points  located  oh  a  boundary  were  obtained  by  approximation  to 
the  boundary  condition  itself.  In  the  present  case,  the  approximation 
to  the  differential  equation  will  be  employed,  and  the  boundary  condition 
incorporated  through  use  of  the  equivalent  symmetry  property. 


Consider  a  typical  point  0  on  the  horizon¬ 
tal  axis  as  shown  in  sketch  (g).  Point  3  Is  a 
fictitious  lattice  point  located  below  the  hori¬ 
zontal  axis  at  9  =  -A,  where  A  is  the  lattice 
interval.  The  finite-difference  approximation  to 
the  differential  equation  (20)  of  the  present 
text  is  given  by  equation  (20)  of  reference  3  as 


/ 

A 


4 • 


♦a  +  *4  -  2Tlo(1ff,i+t,3)-2(l-2Ti0)t'o  =  0  (Al)  j 


where  t|0  is  the  abscissa  of  point  0.  The  symme¬ 
try  property  leading  to  the  boundary  condition  '(38) 


Sketch  (g) 


requires  that  =  t’i,  so  that  for  points  on  the  horizontal  axis  equa¬ 

tion  (Al)  reduces  to 


rz  +  r4  -  -  2(i-2n0)*’0  =  0 


(A2) 
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The  point  at  the  intersection  of  the  horizontal  axis  and  the  shock 
polar  needs  special  consideration.  Sketch  (h)  shows  conditions  at  this 

point.  Here,  as  before,  point  4  is  a  ficti¬ 
tious  point  located  below  the  boundary  sym¬ 
metrical  to  point  3-  It  follows  from  the 
boundary  conditions  (25)  and  (38 )>  both  of 
which  must  be  satisfied  at  the  point  0,  that 
the  first  derivatives  in  the  coordinate 
directions  are  both  zero  at  that  point.  On 
the  basis  of  this  fact,  if  the  function 
T|fi(rj,e)  is  expanded  in  a  two-dimensional 
Taylor's  series  about  point  0,  the  following 
finite-difference  relations  for  the  second 
derivatives  are  easily  obtained: 


A 


Sketch  (h) 


1  \lr  > 

2  0 


a2*' Jo  =  w'i  -r*’i 

A2*’eelo  -  2*’3  -  2*'o  -  ^A2*' Jo 

Here  the  symmetry  property  about  the  horizontal  axis  has  been  used  to 
equate  to  Substitution  of  these  relations  into  equation  (,  Oj 

for  q  =  -1  leads  to  the  following  finite-difference  equation  for  the 

point  0: 

- 1  (1-2 *“)♦•*  +  W'3  -  [4  +  1  d-2ia)  J ♦■0  -  0  <«> 

Finite -Difference  Equations  Special  to  ♦  'p 


The  only  finite-difference  equation  special  to  the  problem  for  t'p 
is  the  one  used  to  terminate  the  field  of  computation  at  some  vertical 
line  on  the  left.  As  in  the  corresponding  work  for  ¥,  this  equation  is 
derived  from  an  asymptotic  solution  of  the  boundary-value  problem  valid 
for  large  negative  values  of  q.  The  derivation  is  parallel  to  that 
described  in  detail  on  pages  l6  and  IT  of  reference  3* 


/r 

\fOl 


0 

for  ?f- 


%(v.o)=° 


0=0 


The  boundary  conditions  which 
must  be  satisfied  by  t'B  at  large 
negative  values  of  q  are  shown  in 
sketch  (i).  The  shaded  section  shows 
the  anticipated  variation  of  t’p 
for  constant  q.  A  solution  of  the 
differential  equation  which  satisfies 
the  given  boundary  conditions  is 


Sketch  (i) 
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^bC^0)  "  X  °n  cos ( 2ew )  x  ‘v^  Kl/3 [ (“2t^)3/2 

n=i 


where  K1/3  is  the  modified  Bessel  function  of  the  second  kind  of 
order  1/3  and  the  Cn  are  appropriate  constants.  If  only  the  leading 
term  of  this  solution  is  used  and  the  Bessel  function  is  replaced  by  the 
first  term  of  its  asymptotic  expansion,  there  results 


As  in  the  earlier  work,  let  A  denote  the  lattice  interval  and  p  some 
large  negative  value  of  q  such  that  A/ j3  « 1 .  It  then  follows  from  the 
foregoing  solution  that,  to  a  first  order  and  for  a  given  value  of  0, 


*'B(-P-A,0) 

+tb(-P>6) 


(a4) 


By  substituting  this  relation  into  equation  (Al),  a  finite-difference 
equation  can  be  obtained  which  is  valid  for  points  on  the  line  q  =  -J3 
and  does  not  include  any  points  to  the  left  of  this  line  (cf.  equation  (22) 
of  reference  3)« 


Finite -Difference  Equations  Special  to 

The  only  equations  special  to  the  problem  for  y arise  as  a  con¬ 
sequence  of  the  condition  along  the  upper  boundary,  where  the  values 
of  are  prescribed  as  a  function  of  q.  Along  most  of  the  boundary, 

this  condition  can  be  met  by  substituting  the .prescribed  values  directly 
into  finite-difference  equations  of  the  type  (Al)  for  points  one  interval 
below  the  boundary.  Because  of  the  nature  of  the  boundary  values 
near  q  =  0,  however,  some  change  from  previous  procedures  is  necessary 
in  the  vicinity  of  the  shoulder.-  Modification  is  also  required  in  the 
equations  used  to  terminate  the  field  on  the  left. 

Points  near  shoulder  of  wedge.-  From  the  known  behavior  of  ?  in 
the  vicinity  of  the  shoulder  (see  Guderley's  results,  reference  13,  for 
the  flow  around  a  convex  corner),  it  can  be  shown  that  the  variation 
of  t'jv  along  the  upper  boundary  near  q  •=  0  must  be  of  the  form 

tTA(n^w)  =  *e(v0w)  =D  (-n)1/2  (a5) 

where  D  is  a  constant  of  proportionality.  '  A  singular  solution  of  the 
differential  equation  (20)  which  is  valid  in  the  vicinity  of  the  shoulder 
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and  which  satisfies  the  boundary  condition  (A5)  is  also  obtainable  from 
Guderley's  results.  This  solution  is,  in  the  present  notation. 


t'As(n^)  =  d(-ti)1/2  (i-C)1/6  f(-  i,  i,  h  JL) 


(A6) 


where  F  is  the  hypergeometric  function  and  £  =  £(q,@)  is  defined  by 


9  (Qw -ey 


£  =  2 

8  q3 

Equation  (A6)  is  suitable  for  use  near  the  upper  boundary  (0  =  0W,  C  =  0). 
Near  the  sonic  line  (-q  =  0,  £  =  -«»)  the  following  alternate  form  is 
available: 

-iVe 


♦a, 


C~T») 


1/2 


(-Ti)3  +  |  (ev-ef 

F(  i,  1 


TP  f  1  1  2.  _1_ 

'  6’  3’  3*  l-t 


22/3(l-01/6  3  3  l-t. 

If  equation  (A7)  is  evaluated  on  the  sonic  line,  there  results 

*as  <°’e>  '  f?75  D  f0*-8)173 


(A7) 


(AB) 


This  result  is  in  agreement  with  equation  (32),  which  was  developed  from 
other  considerations.  It  is  apparent  from  equations  (A5)  and  (A8)  that 
a  solution  for  will  have  a  singularity  in  the  first  derivatives  at 

the  point  q  =  0,  9  =  0W.  .  / 

Because  of  the  foregoing  singularity,  a  direct  numerical  calculation 
of  ty'A  might  be  expected  to  run  into  difficulty  in  the  vicinity  of  the 
shoulder.  Attempts  along  these  lines  lead,  in  fact,  to  the  unlikely 
result  of  negative  lift  over  a  small  region  of  the  profile  just  forward 
of  the  midchord.  Reductions  of  the  lattice  interval  to  quite  small  values 
served  merely  to  decrease  the  extent  of  this  region.  This  is  in  contrast 
to  the  previous  work  for  ^  (and  for  as  well),  in  which  the  singu¬ 

larity  at  the  shoulder  appears  in  the  second  derivatives.  In  that  case, 
a  sufficiently  accurate  solution  for  the  unknown  function  could  be  obtained 
by  direct  calculation.  In  the  present  work,  it  was  found  necessary  to  , 
subtract  out  the  singularity  in  the  first  derivatives  according  to  the 
following  procedure: 


Let  a  function  S^'a  be  defined  such  that 

8*'a  =  t'A-rAs 
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Singular  solution  of  the  type  given  by  equations  (A6) 

r s  s 

be  reiiace4  *>  a  S-  ™ 

is  iinn*  \t>a  it  power  vaflatlon  of  the  form  given  by  equation  (A5).  This 

S  boundary^  %  Tf f"  *?  1“  ' if ' “  ™*  «*>  vithin  tiis 'length 
ary  v  A  v TJ , aw J  -  V^T],^)  or  8^*A  =0.  On  this  basis  a 

boundary-vaiue  /robiem  for  6**A  can  be  defined  for  a  small  region  near 

the  shoulder  as  shown  in  sketch  (j).  The  g 

problem  for  8*r'A  within  this  region  is  - 1 - - - -  B 

solved  jointly  with  the  problem  for  ^'a  £f%'=0  ] 

in  the  remainder  of  the  field.  The  two  '  ! 

regions  are  fitted  together  by  the  use  of  I  r9(X<ir') 

overlapping  lattices,  much  as  is  done  in  I  (S^VkJ /T- 7&Sd4sO  H 

the  case  of  a  graded  mesh  (see  reference  15) .  f  9*  ' 

The  only  difference  is  that  equation  (A9)  - 

must  now  be  utilized  to  make  the  transition  8#=#-*,  1 

between  the  two  lattices  at  all  their  common  Sketch  ( i ) 

points.  It  is  seen  from  sketch  (j)  that  ° 

iSSijTJth  tf'A  °n  b0ththe  uWer  boundary  and  sonic  line  are 
identical  with  the  corresponding  conditions  for  f.  The  finite -different 

SctlTf ~  S'  CalC“latI-  °f  8*’A  =an  therefore^be  tSn  ov^r  ““ 
directly  irom  the  previous  work. 

ora  .As  nearlyn  as  one  can  Judge  from  experience  with  various  lattice 
vSl  23cJli^-  °btainad  by  the  foregoing  process  are  guantitatiyely  as 

S£r“?  0f%T-^g  ss: s&d r- 

errors  from  this  source  are  probably  small.  Physical  plane. 


Sketch  (j) 


Points  far  to  the  left.-  The 
boundary  conditions  for  at 

large  negative  values  of  r\  are 
shown  in  sketch  (k).  From  the 
asymptotic  solution  for  the  basic 
problem  (equation  (21)  of  refer¬ 
ence  3),  the  expression  for  i|r'A 
along  the  upper  boundary  is  found 
to  be 


Sketch  (k) 


ho 
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^’ACn^w)  -  t0(n^w)  =  E  (-ri)"1/4  exp  (-2t])3/2  (AlO) 

where  E  is  a  constant. 

Because  of  the  nature  of  the  boundary  condition  (AlO),  it  is  not 
possible  to  write  an  asymptotic  solution  for  ty*A  for  large  negative  tj 
in  a  single  term.  For  this  reason,  the  procedure  previously  used  to 
terminate  the  field  of  calculation  at  some  location  on  the  left  cannot 
be  applied  in  the  present  case.  An  alternative,  procedure,  somewhat  more 
arbitrary  in  nature,  can  be'  devised  by  writing  tyTA  i*1  form 


*’A(ii,e)  =  rA(n,0w)  +  srA(T],e)  (All) 

% 

where  5t’A  is  now  defined  by  S^’a “  ^’a  (’Ij^w)  (see 
sketch).  The  attenuation  of  ?*A  in  going  from  a  point  at  h  =  ”0 
to  a  point  at  t)  =  -0  -  A  is  then  found  by  treating  each  of  the  terms 
in  equation  (All)  as  an  independent  quantity.  The  attenuation  of 
y  (t\,Qv)  is  found  from  equation  (AlO)  by  a  procedure  similar  to  that  used 

in  obtaining  equation  (A4)  .  The  result  is 

rA(-p-A,ew)  =  (1  -  ^)exP  (-  ^  ^2p).  *fA  ^Mw)  (A12) 

To  obtain  a  corresponding  equation  for  5ty'A  ,  it  is  assumed  that  for  a 
given- value  of  0  this  quantity  attenuates  in  the  same  manner  as  was 
previously  found  for  B'  Cto®  thus  has  from  equation  (A4) 

S**a(-0-A,0)  =  (l  ~  exp  ^  /2^&*’A(-P,0)  (M3) 

Substitution  of  these  expressions  in  equation  (All)  for  .  tj  =  -0 
gives  finally 

rk(-?-A,e)  -  (1-4)  (-  ^•Tw)  x  . 

\\ew  (-^y25)-i]+'A(-P-e«)  +  +-A 

L  L  \  20^  '  J 

Since  ^’a  (-0,0™-)  is  a  known  quantity  for  any  given  value  of  P,  this 
equation  can  be  used  to  terminate  the  field  of  calculation  in  the  same 
manner  as  was  done  with  equation  (A4) .  The  'considerable  element  of  . 
arbitrariness  in  the  derivation  of  equation  (Al4)  can  be  tolerated  since 
the  over-all  solution  is  insensitive  to  changes  in  the  left-hand  portion 
of  the  field. 


-  A 


(Al4) 
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Solution  of  Finite -Difference  Equations 


The  'techniques  used  to  obtain  a  solution  of  the  finite -difference 
equations  for  \j and  ty'g  were  the  same  as  those  described  in  refer¬ 
ence  3  for_the  basic  solution  of  f.  In  general,  the  graded  lattice  as 
used  for  if  (see  fig.  2  of  reference  3)  was  suitable  for  the  solution 
of  t’jj.  For  however,  different  gradations  were  necessary  with 

the  smallest  lattice  spacing  being  used  near  the  shoulder  (point  B  of 
sketch  (e)).  The  value  of  at  the  intersection  of  the  shock  polar 

and  the  sonic  line  was  chosen  as  10,000  so  that  the  previously  obtained 
values  of  if  could  be  used  to  provide  the  initial  guess  for  ty'p. 

In  the  course  of  the  present  work,  a  useful  technique  was  found  for 
locating  regions  of  relatively  large  error  in  the  numerical  solution. 

By  use  of  one  form  of  Green’s  theorem  plus  the  differential  equation  (20 )s, 
it  can  be  shown  that  around  any  contour  enclosing  a  region  in  which 
equation  (20)  is  satisfied  the  following  relation  must  hold: 

.  .  j»  ■ 

f  (2qf*edri  +  V'^de)  =  0  (A15) 

In  a  numerical  solution  the  line  integral  in  equation  (A15)  will  not, 
except  by  rare  coincidence,  be  precisely  zero  around  any  given  contour. 

The  amount  by  which  it  differs  from  zero  may  be  taken  as  a  rough  measure 
of  the  adequacy  of  the  numerical  solution  over  the  region  within  the 
contour.  If  the  entire  field  of  calculation  is  subdivided  into  a  number 
of  contiguous  regions,  it  is  thus  possible,  by  evaluating  the  integral 
around  each  of  the  enclosing  contours,  to  locate  regions  within  which 
the  error  is  relatively  high.  The  solution  in  these  regions  can  then  be 
improved  by  advancing  locally  to  a  finer  mesh.  This  technique  was  found 
to  be  of  great  help  in  the  present  work.  It  would  probably  be  useful  in 
other  elliptic  boundary-value  problems  for  which  a  relation  analogous  to 
equation  (A15)  can  be  obtained. 


NACA  W  2832 


4-2 


APPENDIX  B 

CALCULATION  OF  FLOW  OVER  REAR  WEDGE 
IN  PHYSICAL  PLANE 


The  procedure  used  to  calculate  the  flow  over  the  rear  wedge  has 
been  outlined  in  the  section  METHOD  OF  SOLUTION.  The  fundamental  opera¬ 
tion  is  to  determine,  by  stepwise  methods,  the  initial  rate  of  movement 
of  the  known  intersection  points  in  the  basic  characteristics  net.  The 
methods  which  are  used  depend  on  the  fact  that  these  points  are,  by 
virtue  of  the  basic  characteristics  construction,  points  of  fixed  r j,0 
(cf.  equations  (55)  and  (57)  of  reference  3).  , 

The  first  step  is  to  determine  the  initial  rate  of  movement  of 
those  points  at  which  the  Mach  lines  of  the  basic  characteristics  net 
meet  on  the  sonic  line.  For  this  purpose,  consider  equations  (48),  - 

which  give  the  initial  rate  of  movement  of  a  general  point  of  fixed  rj,e. 
If  these  equations  are  specialized  to  apply  to  points  on  the  sonic  line, 
the  following  relations  are  obtained: 

(Bla) 


1  P 9 

X'(0 ,0)  =  -=“  / 


Y’(o,e) 


&v?/S 


r 


(Bib) 


To  write  equation  (Bla)  the  path  of  integration  in  equation  (4-8a)  is 
taken  along  the  upper  boundary  from  0  to  B  (see  sketch  (e)  on  page  21) 
and  thence  downward  along  the  sonic  line.  The  contribution  of  the 
portion  from  0  to  B  is  zero  by  virtue  of  the  condition  (37)*  In  apply¬ 
ing  these  equations,  the  value  of  Iw  is  known  from  the  basic  solution. 
The  integral  in  equation  (Bla)  is  evaluated  by  mechanical  integration 
of  a  curve  of  numerically  determined  derivatives.  Proper  allowance  is 
made  for  the  singularity  at  the  shoulder  by  integrating  the  singular 
solution  analytically.  The  component  rates  of  movement  of  the  sonic 
point  at  the  shoulder  are  both  seen  to  be  zero-. 


The  next  step  in  the  solution  is  to  calculate  the  rate  of  movement 
of  intersection  points  downstream  of  the  sonic  line.  This  is  done  by 
proceeding  stepwise  along  consecutive  downgoing  characteristics. 
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Consider  three  typical  net  points  as  shown  in  sketch  (l)  (cf.  also 
fig.  4  of  reference  3)-  The  dashed  lines  represent  the  original  posi¬ 
tion  of  the  Mach  lines  through  points  a,  b,  and  c,  and  the  solid  lines 
represent  their  displaced  posi¬ 
tions  corresponding  to  a  small, 
finite  angle  of  attack  a. 

Since  the  intersection  points  in 
the  Mach  net  are  points  of  fixed 
r 1,0,  the  components  of  their 
displacement  are  given  by  aX' 
and  aY* .  The  slope  of  each  segme 
of  Mach  line  is  taken,  in  accord 
with  the  procedures  of  refer¬ 
ence  3,  as  the  average  of  the 
slopes  calculated  at  the  two  end 
points.  The  slope  calculated  at 
each  end  point  depends,  in  turn, 
only  on  the  value  of  rj  at  that 
point  (cf.  equation  (54)  of 
reference  3) • 

* 

It  is  desired  now  to  deter¬ 
mine  X*  and  Y*  at  point  c  in  Sketch  (l) 

terms  of  X'  and  Yr  at  points  a 

and  b.  Since  the  value  of  T)  at  a  given  net  point  is  the  same  in  the 
displaced  and  undisplaced  positions,  it  follows  from  what  has  been  said 
above  that  each  Segment  of  Mach  line  must  retain  its  original  slope 
after  displacement.  If  this  slope  is  denoted  by  m,  the  following 
relations  are  then  readily  obtained: 


.  Y’a -Y*b+mjjCX*b -fiaQX'a 

x  c  =  - ; - - - -  (B2a) 

®bc  "  ®ac 


^bc^'a~^acT,b+^ac?bc(^,b  ~X*a) 
™bc  “  ™ac 


(B2b) 


With  these  relations,  it  is  a  simple  matter  to  calculate  the  initial 
rates  of  movement  of  successive  net  points  on  consecutive  downgoing 
characteristics.  For  the  first  characteristic  to  be  considered,  point 
b  is  taken  at  the  shoulder  of  the  profile,  where  X'  and  Y»  are  both 

zero.  Thus,  X'c  and  Y’c  for  net  points  on  this  characteristic  can  be 

determined  solely  in  terms  of  X’a  and  Y»a  and  the  slopes  maC  and  mbc* 

For  the  remainder  of  the  downgoing  characteristics,  X’b  and  Y'b  are 
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known  from  calculations  along  the  characteristic  immediately  preceding. 
The  actual  calculations  can  he  carried  out  in  straightforward  tabular 

form. 

The  foregoing  procedure  enables  the  calculation  of  X’  and  Y*  for 
all  net  points  except  the  ones  originally  at  the  surface  of  the  rear 
wedge.  For  these  points,  consideration  must  be  given  to  the  required 
boundary  condition  at  the  surface.  This  boundary  condition  is 

0(X,+O;a)  =  -  (Gy+a)  (®3) 

from  which  it  follows  that 

0»(X,+O)  =  -  1  (B4) 

The  problem  now  is  to  determine  X*  and  Y*  at  the  surface  of  the 
wedge  in  such  a  way  that  equation  (b4)  is  satisfied.  To  do  this  equa¬ 
tion  (53b)  is  first  specialized  to  the  surface  of  the  wedge,  where  it 
is  readily  shown  that  X0  «  Ytj  =  0.  In  view  of  condition  (b4),  there 
results  * 

Y*(q,-0w)  =  Y6(ti,-0w)  (B5) 

The  value  of  Y*  at  points  originally  on  the  surface  of  the  wedge  is 
thus  fixed  directly  by  the  basic  solution.  The  corresponding  value  of 
X'  can  be  found  from  a  construction  analogous  to  that  of  sketch  (l) 
and  is 

x,c  =  Y'a-n  -5acX'a  (b6) 

-5ac 


The  point  c  is  now  the  point  originally  on  the  surface  of  the  wedge 
(i.e.,  Y'c  is  as  given  by  equation  (B5))>  and  the  remaining  notation 
is  the  same  as  in  sketch  (l) . 


Application  of  equation  (B5)  requires  the  knowledge  of  Y0(t),-0v), 
which  in  the  case  of  the  wedge  profile  is  equal  to  l/0y.  Evaluation 
of  the  latter  derivative  can  be  carried  out  directly  from  the  basic 
Mach  net,  but  the  procedures  are  cumbersome  and  inaccurate.  A  better 
method  is  to  use  the  equations  of ‘motion  (cf.  equation  (6)  of  refer¬ 
ence  4)  to  express  0y  in  terms  of  fj^r.  Following  this  procedure,  one 
obtains  finally 


Yq ( 4 , ”®w) 


2fj(X,0)fix(X,+0) 


(B7) 
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The  quantities  f]  and  which  appear  here  are  easily  evaluated  from 
the  basic  solution  for  the  chordwise  distribution  of  fj. 

The  preceding  equations  enable  the  calculation  of  the  initial  rate 
of  movement  X*  for  points  originally  on  the  surface  of  the  rear  wedge. 
The  final  step  is  to  determine  the  corresponding  distribution  of  lift.. 
For  this  purpose,  equation  (53a)  is  specialized  to  points  on  the  rear 
wedge  to  obtain 


y(x,40)  =  -  x'(n,-ew)/x^(fi,-ew) 

which,  in  view  of  the  boundary  conditions,  can  be  shown  to  be  equivalent 
to 


V(x,+o)  .  -  x»(r],-ew)fjx(x,+o)  (b8) 

The  distribution  of  lift  is  then  obtained  from  equation  (57) • 
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APPENDIX  C 

SOLUTION  OF  PROBLEM  FOR  COMPLETELY  SUPERSONIC  FLOW 

Calculation  of  Lift -Curve  Slope  and 
Center  of  Lift 


If  conditions  are  such  that  21/3  =  1.260  (corresponding 
to  ©wS •  equation  (13)),  then  the  basic  flow  over  the  profile  at 
zero  angle  of  attack  is  completely  supersonic.  The  solution  for  the 
lift -curve  slope  and  center  of  lift  at  a  vanishingly  small  angle  of 
attack  can  then  be  carried  out  analytically  as  follows: 

Consider  a  completely  supersonic  flow  about  the  double-wedge  profile 
at  a  small  angle  of  attack.  In  the  physical  plane  the  flow  field  has 
the  well-known  appearance  shown  on  the  left  in  sketch  (m) .  The  corre¬ 
sponding  hodograph  of  the  flow  along  the  upper  surface,  in  terms  of  the 

ff 


normalized  small-disturbance  variables  T]  and  0,  is  shown  on  the  right. 
The  quantities  0W  and  a  are,  as  before,  the  half -angle  of  the  wedge 
and  the  angle  of  attack  (also  normalized) .  Except  for  a  small  range 
of  just  above  1.260  (see  below),  flow  conditions  must  be  constant 

along  each  of  the  segments  1  and  2  of  the  upper  surface.  In  the  hodo¬ 
graph  each  of  these  segments  is  thus  represented  by  a  single  point 
located  as  shown.  It  is  apparent  that  for  a  given  value  of  0V,  the 
speeds  r]  and  r)2,  which  are  the  primary  unknowns  in  the  problem,  are 
functions  solely  of  the  angle  of  attack  a. 
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To  find  the  lift-curve  slope  and  center  of  lift  it  is  necessary 
first  to  find  the  derivatives  Vi  =  (dT]1/da)a_0  an<^  7  *2  =  (dT]2/ <ia) a=o * 
This  can  he  done  with  the  aid  of  the  equations  for  the  shock  polar 


Compare  equation  (53)  of  reference  3- 
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To  find  t]’2,  equation  (C2)  for  the  downgoing  characteristic  is 
first  specialized  so  as  to  pass  through  the  point  1.  This  gives 

0  >  (e„-a)  +  ^  (%3/2  -  13/2) 

Substitution  of  the  boundary  condition  92  =  -  9V  -  a  then  provides 
the  result  that 


2/3 


Taking  the  derivative  with  respect  to  a,  one  obtains  finally  at  a  =  0 


-  1/2 


C 


,3/2  +  J_  ^/3 

J2- 


*1  1”/”+  ~T=-ev) 


Vi 


(c  6) 


where  T]^  is  given  by  equation  (c4)  and  t|  by  equation  (C5). 

Since  the  value  of  tj'  is  constant  on.  each  segment  of  the  profile, 
the  lift-curve  slope  is  easily  found  from  equation  (57)  and  is 

[(7+i)(t/c)]1/3  f-^f)  =  2(2ew)1/3(v1+T]»2) 

x  ■  a=o 


Substitution  from  equation  (C6)  gives 

«7+l)(t/c)J1/3  (^)  -  2(2ew)l/Vi 

x  /a=o 


T?  X/2 


1  + 


5. 3^4:0^ 


(V 


(07) 


The  moment-curve  slope,  for  moments  taken  about  the  leading  edge,  is 
found  to  be 

[(7+l)(t/=)]1/3(^)  -  (20»f/3U,1+  3i'2) 

x  'a=o 


or 


|V3/±m\ 


((r.ixt/cr^^  -  -  I 


1  + 


3V1/2 


nx3/s  + 


®wy 


\  1/3 


(os) 
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The  position  of  the  center  of  lift  is  given  accordingly  by 


x  +  3f^i1/g 

( dcm/ da)  2.  [^3/2+  (3 /*/2)6v]1/3 

(dcj/da)a_0  ^  i  +  %1/g _ 

[r\!3/£+  (3/i)^] 1/3 


In  equations  (C7)  and  (C8),  the  first  term  inside  the  brackets  repre¬ 
sents  the  contribution  of  the  front  wedge ,  the  second  term  that  of  the 
rear. 


Equations  (C7)  and  (C9)  are  the  basis  for  the  curves  shown  in 
figures  5  and  7  for  values  of  £o^1.260.  The  results  show  certain 
curious  features  when  the  flow  over  the  front  wedge  is  just  sonic 
(i.e.,  ^=0,  0^=1,  6o=1.260).  These  are  as  follows: 

(a)  The  lift  contributed  by  the  rear  wedge  is  zero  (see  equa¬ 
tion  (C7). 

(b)  The  center  of  lift  is  at  the  quarter-chord  point  (follows 
from  statement  (a)  plus  the  condition  of  uniform  lift  on  the  front 
wedge;  see  also  equation  (C9)). 

(c)  The  rate  of  change  with  respect  to  |  is  infinite  both  for 
the  lift-curve  slope  of  the  complete  profile  ana  for  the  position  of 
the  center  of  lift  (follows  from  differentiation  of  equations  (C7) 
and  (C9)). 

These' results  are  associated  in  every  case  with  the  behavior  of  the  lift 
calculated  for  the  rear  wedge. 


Estimation  of  Lower  Limit  for  Constant  Speed 
Along  Rear  Wedge 


The  features  just  enumerated,  though  having  a  certain  curiosity  in 
themselves,  cannot  be  accepted  as  completely  correct.  Because  of  inter¬ 
action  effects  between  the  shock  wave  from  the  bow  and  the  expansion 
fan  from  the  shoulder,  the  fundamental  condition  of  constant  speed  at 
the  surface  of  the  profile  will  not  be  satisfied  along  the  rear  wedge 
•until  the  value  of  |0  is  somewhat  greater  than  1.260.  Until  then, 
disturbances  reflected  from  the  shock  wave  will  reach  the  rear  wedge 
and  cause  a  slight  decrease  in  speed  toward  the  trailing  edge.  This 
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effect  will  cease  when  the  forwardmost  reflected  Mach  wave  just  touches 
the  trailing  edge.  The  exact  value  of  |Q  at  which  this  condition  will 
he  met  is  difficult  to  determine.  An  upper  hound  can,  however,  he  esti¬ 
mated  as  follows: 

Consider  the  "basic  flow  field  (o=0)  over  the  upper  half  of  the 
profile  when  the  first  reflected  Mach  wave  just  strikes  the  trailing 
edge.  Sketch  (n)  shows  such  a  flow  field  as  it  would  appear  in  tran¬ 
sonic  similarity  form  (cf.  pp.  12-13  of  reference  2).  In  drawing  the 
sketch  a  special  assumption  has  been  introduced  beyond  those  implicit 
in  the  small-disturbance  theory;  namely,  that  the  first  reflected  Mach 


Sketch  (n) 

wave  is  straight  and  has  an  angle  of  inclination  p  equal  to  that  of 
the  first  wave  in  the  expansion  fan.  With  this  assumption,  the  corre¬ 
sponding  value  of  !0  is  ‘easily  determined.  Since  the  reflected  wave 
must  actually  he  curved  downstream,  the  value  so  determined  will  he 
greater  than  the  correct  value  for  the  required  condition. 

On  the  basis  of  sketch  (n),  the  following  equation  can  he  written 
between  the  shock  angle  X  and  the  Mach  angle  p: 

tan  X  =  i  tan  p  (CIO) 

A  relation  between  the  shock  angle  X  and  the  speed  fj1  in  the  region 
behind  the  shock  can  be  obtained  from  equation  (C3)  and  the  known  prop¬ 
erties  of  the  shock  polar.  The  result  is 
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tan  \ 


(2MV3 

'</l+f)1 


An  analogous  expression  for  the  Mach  angle  q  is  given  by  equation  (5^-) 
of  reference  3  and  is 


tan  |i 


(^)1/3 

*/  2^1 


Substitution  of  these  relations  into  equation  (CIO)  and  solution  for  fL 
gives 


The  accompanying  value  of  GWf  found  from  equation  (C3)  with  a=0,  is 


ev  =  0.9685 

This  corresponds,  according  to  equation  (13),  to 

£0=  1-287  (Cll) 


Thus,  for  values  of  between  1.260  and  some  limit  less  than  1.287, 

the  results  of  equations  (C7),  (C8),  and  (C9)  are  not  exact  insofar  as 
the  contribution  of  the  rear  wedge  is  concerned.  It  can  be  reasoned 
that  in  this  range  an  exact  solution  would  indicate  more  lift  for  the 
rear  wedge  than  does  the  present  analysis. 
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Figure  4.-Chordwise  distribution  of  lift. 
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Figure  5.-  Generalized  lift-curve  slope  as  a  function  of  transonic  similarity  parameter. 
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□  Present  numerical  analysis  -  Transonic  sma/J  disturbance  theory 

o  Guderley  and  Yoshihara  (Ref.  9)  —  —  Linear  theory 
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Figure  7.  -  Center  of  lift  as  a  function  of  transonic  similarity  parameter. 
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